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ANALYTICITY OF HOMOGENIZED COEFFICIENTS UNDER BERNOULLI 
PERTURBATIONS AND THE CLAUSIUS-MOSSOTTI FORMULAS 

MITIA DUERINCKX AND ANTOINE GLORIA 


This paper is concerned with the behavior of the homogenized coefhcients associated with some 
random stationary ergodic medium under a Bernoulli perturbation. Introducing a new family of 
energy estimates that combine probability and physical spaces, we prove the analyticity of the per¬ 
turbed homogenized coefficients with respect to the Bernoulli parameter. Our approach holds under 
the minimal assumptions of stationarity and ergodicity, both in the scalar and vector cases, and gives 
analytical formulas for each derivative that essentially coincide with the so-called cluster expansion 
used by physicists. In particular, the first term yields the celebrated (electric and elastic) Clausius- 
Mossotti formulas for isotropic spherical random inclusions in an isotropic reference medium. This 
work constitutes the first general proof of these formulas in the case of random inclusions. 
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1. Introduction 

Let p = {qn)n be a stationary and ergodic random point process in of intensity unity. To each 
point Qn we associate a family of independent Bernoulli variables [0,1] D p i—)• 6^^ that takes value 1 
with probability p and value 0 with probability 1 — p. Given a, /3 > 0, we define a family of random 
matrix fields on as follows: 

CX) 

A^P\x) := ald + '^U;P\l3 - 

n=l 

where B[qn) denotes the unit volume ball centered at q^ (above we have assumed that the balls 
B{qn) are disjoint for simplicity of the discussion). The random matrix fields A^'P^ are stationary and 
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ergodic, so that the standard stochastic homogenization theory [27] yields the existence of associated 
homogenized matrices For p small, we expect to be a perturbation of order p of the 

unperturbed medium aid. Indeed, in the late 19th century, Clausius (1879) and Mossotti (1850) 
independently proposed the following expansion in the scalar case (also known as the Maxwell relation 
[24]): 


A 


(p) _ 


horn 


= ald+' 


ad(/3 — a) 


■ldp + o(p). 


( 1 . 1 ) 


/3 + a{d — 1) 

The rigorous proof of this statement and of its counterpart for linear elasticity has remained a chal¬ 
lenge since then. The first justification of the (electric) Clausius-Mossotti relation is due to Almog in 
dimension d = 3, whose results in [1, 2], combined with elementary homogenization theory, precisely 
yield (1.1) (the convergence rate obtained in [2, Theorem 1] is lost when combined with homogeniza¬ 
tion). The proof is based on (scalar) potential theory and crucially relies on the facts that d = 3, that 
Aip) is everywhere a multiple of the identity, and that the inclusions are supposed to be spherical 
and disjoint. Another contribution is due to Mourrat [25], who considered a discrete scalar elliptic 
equation instead of a continuum elliptic equation for all d > 2. In the case treated in [25], A^P'^ is a 
discrete set of i.i.d. conductivities (5 and a with probabilities p and 1—p. The extension of Mourrat’s 
results to the present continuum setting (which is made possible by the recent contributions [ 12 , 16]) 
would yield the improvement of ( 1 . 1 ) to 

ad[l3 — a) 


\{p) 

^hom 


= ald-|- 


/3 -|- a{d — 1) 


I6.p + 0{p '’') 


( 1 . 2 ) 


for any 7 > 0. The assumptions of [25] are however very stringent, and typically cover the case of 
a hardcore random Poisson process (i.e. a Poisson point process where overlapping inclusions are 
deleted). Indeed, Mourrat’s proof crucially relies on quantitative estimates on the corrector obtained 
in a series of works [13, 14, 15, 23] by the second author, Marahrens, Neukamm, and Otto. These 
estimates typically hold under a spectral gap estimate assumption, whence the restriction on the 
random point set. There is a small gap in the proof of (the discrete counterpart of) (1.2) (a quantitative 
estimate of {aKp) — a^(0)[, see [25, (8.2)], is missing to complete the proof of [25, Theorem 11.3]), 
which we think can however be fixed using the quantitative results of [14, 15, 23]. Mourrat also made 
the nice observation that the reference medium needs not be the unperturbed background medium 
(of conductivity a), and essentially proved that p 1 —)• is on the whole interval [ 0 , 1 ] for 

all 7 > 0 in the discrete setting. 


The idea of perturbing a non-uniform reference medium first appeared in the work [4] by Anan- 
tharaman and Le Bris, who considered the perturbation of a periodic array of inclusions by i.i.d. 
Bernoulli variables (the inclusions are independently deleted with probability p). The corresponding 
matrix field A^P^ is then a random ergodic field with discrete stationarity. As above, for p small, one 
expects to be a perturbation of (the homogenized coefficients of the unperturbed periodic 
medium) of order p. Anantharaman and Le Bris considered the approximation of obtained 
by periodizing the random medium A^^) on a cube of size L. The qualitative homogenization theory 
ensures that almost surely the random approximation A^^ of A^^^ converges to (note that 

A^^ is deterministic and coincides with A^^^ for all L G N). Although not formulated this way, they 
essentially proved that the approximation p 1 —)• E[A^^] is at p = 0 , and obtained bounds on the 
first two derivatives 9pE[A^^]|p=o '9p®'[A^^]|p=o that are uniform in L. This is however not quite 

enough to prove that p 1 —)■ is itself (or C^) at zero. 

In the present contribution we shall prove in a very general setting (which includes both the 
examples studied by Mourrat and by Anantharaman and Le Bris) that the map p 1 —)• is analytic 

on [0,1] (see Theorem 2.1). Our result holds under the mildest statistical assumptions on the reference 
medium x e->• A(^)( 3 :) and on the point process, that is (discrete or continuum) stationarity and 
ergodicity. We also make the (crucial) technical assumption that the number of intersections is 
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uniformly bounded (see however Remark 2.7 for the specific example of a Poisson point process). We 
also believe that a suitable adaptation of our arguments may allow to treat the case when the Bernoulli 
law is replaced by more general laws as considered in [5, Section 3]. Although our results are much 
stronger than those of Mourrat [25], our proof was mainly inspired by the ingenious computations of 
Anantharaman and Le Bris (see in particular [3, Proposition 3.4]), and only relies on soft arguments. 
In particular the crucial ingredient of our proof is a new family of energy estimates that combine 
both physical and probability spaces, see Proposition 4.6 below. Since the proof only uses ingredients 
that are available for systems, our results hold not only for scalar equations, but also for uniformly 
elliptic systems and for linear elasticity. In the case of an isotropic constant background medium 
perturbed by randomly distributed isotropic spherical inclusions, this proves the celebrated (electric) 
Clausius-Mossotti formula with an optimal error estimate (see Corollary 2.4), 


A 


ip) 

horn 


a Id 


adiji — a) 
(3 + a{d — 1) 


Idp + O(p^), 


(1.3) 


as well as its elastic counterpart (see Corollary 2.5), under the weakest assumptions possible. 

Our proof makes use of the standard modification of the corrector equation by a massive term of 
magnitude and the derivatives of with respect to p are given by the limits as T f oo of the 

derivatives of a deterministic approximation of {VT plays a similar role as the period L in 

[4]). Interestingly, in the case when the inclusions are disjoint, the fact that p i—)• is on [0,1] 
can be obtained as a corollary of the (classical) energy estimates of [4, 3] (applied to A^^ instead of 
A^^), further using that they hold for any p S [0,1] by Mourrat’s observation, see Section 3. The 
proof that p i—)• a!^}^ is analytic is however more subtle and is based on the new family of energy 
estimates. 

There are two motivations to go beyond First, this gives a definite answer to the maximal 

iv) 

regularity of the map p i—)• A^^^. Analyticity was indeed conjectured in applied mechanics and applied 
physics, see for instance [32, Chapters 18 & 19]. This analyticity result contrasts very much with the 
corresponding regularity of a similar model originally studied by Maxwell [24] and Rayleigh [29], see 
also [19, Section 1.7] and [6]. The latter consists of a homogeneous medium periodically perturbed 
by inclusions of volume p located at integer points of (note that the perturbed medium is still 
Z'^-periodic), in which case the associated map p i—)■ is in dimension d, but not more. 

The second motivation stems from numerics. Indeed the main motivation for [4, 5] is to exploit the 
perturbative character of for p small (seen as a model for defects) and use the first two terms of 

a Taylor-expansion at zero as a good approximation for , the accuracy of which can be optimally 
quantified by Theorem 2.1 below. Note that Legoll and Minvielle [22] made an original use of this 
approximation of in a control variate method to reduce the variance for the approximation of 

homogenized coefficients. 

As emphasized in [4], a second natural question is to quantify the convergence speed of these 
approximations of the derivatives of (which are indeed computable in practice). As opposed to 

Theorem 2.1, which is a qualitative result, establishing such a convergence speed requires to quantify 
the speed of convergence in the ergodic theorem for stochastic homogenization, which is a quantitative 
result. We give such a result under the assumption that the speed of convergence of At to Ahom can 
be quantified, which typically (but by far not exclusively) holds in the case that the equation is 
scalar and that A^^^ satisfies a spectral gap inequality (as assumed in ]25]) using results of [16] (see 
Corollary 2.6). 


Our approach to prove Theorem 2.1 is constructive and explicit bounds are obtained on the deriva¬ 
tives. We do not know whether there is an abstract alternative to prove analyticity. In [8], Cohen, 
Devore and Schwab obtained a result of the same flavor using a complexification method: they proved 
the analyticity of the solution of linear elliptic PDFs with respect to parameters in the coefficients in 
the framework of a chaos expansion. Their setting is however very much different from the setting of 
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Theorem 2.1. Indeed, as emphasized in Remark 2.8 below, the solution of interest here (the corrector) 
is not even differentiable with respect to the Bernoulli parameter (the homogenized coefficients are 
analytic because of subtle cancellations). 

It is also not clear to us how much the Bernoulli law can be relaxed to allow for correlations. In 
particular it would be interesting to determine whether the regularity of the homogenized coefficients 
depends on the decay of correlations of the generalization of the Bernoulli law. 

For the clarity of the exposition, although our proof of Theorem 2.1 holds in the case of uniformly 
elliptic systems and of linear elasticity (provided the elasticity tensor is uniformly very strongly elliptic, 
as standard in homogenization), for non-symmetric coefficients, and for discrete elliptic equations, 
we use continuum scalar notation and assume the coefficients are symmetric. For non-symmetric 
coefficients, it is indeed enough to consider, in addition to the primal corrector equation, the dual 
corrector equation (associated with the pointwise transpose coefficients), which would only make 
notation heavier. In addition we assume that the coefficients enjoy continuum stationarity (in the 
case of Z'^-stationarity, the expectation would simply be replaced everywhere by the expectation of 
the integral over the unit cube). Note that our result also covers the case of laminates, or more 
generally the case when the heterogeneous coefficients are random in some direction(s) and invariant 
along the other direction(s) (cf. the example of cylindrical fibers considered in [4] and encountered in 
practice). 

The rest of the article is organized as follows. In Section 2 we introduce the main notation and 
state the main results of the paper: the analyticity of the homogenized coefficients with respect to 
the Bernoulli parameter and the validity of the Clausius-Mossotti formulas. In Section 3, we present 
the general strategy of the proof under the additional simplifying assumption that the inclusions are 
disjoint. Section 4 is dedicated to the introduction and proofs of auxiliary results, and in particular 
of the improved energy estimates. The main results are proved in Section 5. 

2. Main results 

2.1. Assumptions. Let A be a random held. We choose a point process p = {qn)n, and random 
bounded inclusions {Jn)n centered at the points To the inclusions {Jn)n we attach i.i.d. Bernoulli 
variables {hn'^)n with parameter p £ [0,1], and we perturb A on if = 1. The only assumptions 
we need here are stationarity and ergodicity, as well as some deterministic bound on the degree of 
intersections between the inclusions. More precise definitions are given below. 

In the rest of this paper, we denote by Br{x) the ball of radius r centered at x in and we simply 
write Br := Br{0) for balls centered at the origin, and set B := i?i(0). We also denote by Q the unit 
cube [—1, \Y. 

Point process. Let p be a (locally finite) ergodic stationary point process on and choose 
for convenience a measurable enumeration p = {qn)'Y=i- open set D of we denote by 

p{D) := i^{qn & D,n £ N} the number of points of p in D. 

Inclusions centered at the point process. Let i? > 0 be fixed. For all n, let J° be random 
Borel subsets J° C Bji{c. M'^) (maybe depending on p = {qn)n)- This defines random bounded Borel 
inclusions := qn + JY We assume that this inclusion process is stationary, in the sense that the 
random set IJ^ Jn is stationary. Moreover, we further assume that the intersections between the 
inclusions J„’s are of degree bounded by some deterministic constant F £ N; by stationarity, this just 
means 

#{n £ N : 0 £ Jn} < F, almost surely. (2T) 

In physics (see e.g. [32, Section 3.1]), the constant F in (2.1) is called the impenetrability parameter: 
different inclusions may penetrate each other but only with the fixed finite maximum degree F. 

As Jn C Bji{qn), assumption (2.1) is trivially satisfied if we assume p{Q) < 9q a.s. (thus forbidding 
arbitrary large clusters in the point process), but it is important to note that the only problem is 
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the possibility of intersections of arbitrary large degree, which has a priori nothing to do with the 
point process p itself. In the case of inclusions with inner radius bounded from below, however, 
assumption (2.1) is equivalent to an assumption on p of the form p{Q) < 9q almost surely. 


Reference random fields. Given 0 < A < 1, denote by M\ the space of uniformly elliptic 
symmetric d x d-matrices M satisfying ^ for all ^ G Let A, A' be two M.\- 

valued ergodic stationary random fields on Note that A and A' do not need to be independent of 
the point process p, and we simply assume that this dependence is local, in the sense that A(0) and 
^'(0) only depend on p via the restriction p\b^ for some given deterministic r > 0. 

Bernoulli perturbation of A. For any fixed p G [0,1], we choose a sequence {bn^)n of i.i.d. 
Bernoulli random variables with = 1] = p, independent of all previous random elements. We 

can now consider the following p-perturbed random field, which is a perturbation of the random field 
A on the inclusions for which bii' = 1: 


= A1 




=e(p) 




where we have set := {n G N : bn'^ = 1}. 

In the case when the inclusions are disjoint, the p-perturbed random field can be rewritten as 
follows: 

•4“ = E 

n 

Moreover, in that case, as A' is allowed to depend (locally) on the inclusion process, the following 
interesting particular example can be considered: choose a sequence (M(^)n of i.i.d. A^;v-valued random 
fields, and define 


A' :- IdlRd\|j^ Jn + X] 

n 

so that takes the form 

■4'’” = E ('’l'”'^'-. + (1 - '>1'”)-^) 

n 


( 2 . 2 ) 


Probability space and product structure. Let us now briefly comment on the underlying 
probability space. Let (fli, Ju,Pi) be a probability space on which the (stationary) random elements 
/O, {Jn)m and A' are deflned. For all p G [0, 1] and re G N, let := G {0,1}}, endowed with 
the trivial cr-algebra , and let be the Bernoulli measure of parameter p on ■ The probabil¬ 
ity space we consider in this article is the product space (12, P) of (12i, J^i,Pi) and 
for all p G [0, 1] and re G N (with the cylindrical cr-algebra). With E, Ei, E^^^ the expectations with 

respect to the measures P, Pi, P^^^, respectively, we have by definition E = Ei Hpefo i] OneN^^^n- 
The independence of the Bernoulli variables, at the origin of this product structure, then takes 
the form: for any integrable random variables Xn^ Vn^ defined on fli x Um,m^n^2}n ^nd 12^^^, 
respectively, 

E[x(f)r?(P)]=E[x(f)]E[#]. (2.3) 

Note that for all p G [0,1] the random field A^'> is defined on 12i x OnsN stationary 

and ergodic for the measure Pi (8) ( 2 )neN^ 2 ^n- such, it can be viewed as a stationary and ergodic 
random field on (12, P). 

Typical examples. One typical example is that of spherical inclusions = B/j(g„) centered at 
the points of any ergodic stationary random point process p with minimal distance bounded away from 
0. In that case, p can be chosen as the hardcore Poisson point process (that is, a modified Poisson 
process for which points that are at a distance less than 2Rq are deleted; see also the hardcore 
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construction in Step 1 of the proof of Theorem 2.1 in Section 5.2) or the random parking measure 
(see [28], and see also [17] for the ergodicity). Instead of spherical inclusions, we can consider more 
general (and random) shapes Jn = <ln + Jni where the J°’s are i.i.d. copies of some random Borel set, 
with J° C Br a.s. 

Another interesting example is when p is a Poisson process (or any other ergodic stationary point 
process) and when is the Voronoi cell at intersected with the ball at of radius R, say. We 
could alternatively choose for the largest ball of radius less than R centered at Qn and completely 
included in the Voronoi cell at In this case, the number of inclusions per unit volume is not 
necessarily uniformly bounded. 


2.2. Notation. We start with the definition of homogenized coefficients, correctors, and approximate 
correctors. We then introduce the crucial notion of difference operators, and conclude with the 
introduction of a notational system for perturbed coefficients which will turn out to be very convenient 
when it comes to the (numerous) combinatorial arguments involved in the proofs. 

Correctors, approximate correctors, and homogenized coefficients. For any (possibly 
infinite) subset FI C N, we define := A + C®, where ;= (A' — A)1lje and := UneE 

these terms, with := {n S N : bn'^ = 1}, we have = A + C^^^\ and we use the short-hand 

notation := and A^P'> := A^^^\ 

For all T > 0, we define the approximate correctors and in direction |^| = 1, associated 
with any field A and any A^, respectively, as the unique solutions in the space {v € : 

sup^-|- |Vup) < oo} of the equations 

^(/)t, 5 -V-A(V</-t,c + 0 = 0, and _ V • + 0 = 0. (2.4) 

These solutions satisfy the following energy estimates (see [16, Lemma 2.7]): 

sup/ (r-VT,cP + |V</>T,d') < 1 , sup/ (T-Vf,^|" + |V0f,/) < 1 . (2.5) 

To shorten notation, we write for For all p E [0,1], as the random field A^P'^ is ergodic 

and stationary, we have (combine for instance the ergodic theorem with [27] in the symmetric case, 
and [11, Theorem 1] in the non-symmetric case) 

limE[[V(/>^^g-V(/>^^^]2] = 0, (2.6) 


(p) 

where Vcjy^ is the gradient of the corrector, i.e. the gradient of the unique measurable random map 
G L^qj,(M'^) solution of the equation 

-V-+0 =0 

on that satisfies = 0 almost surely and such that is stationary and has bounded 

second moment. (Note that exists for any matrix field A if T > 0, whereas (j)^ only exists almost 
surely for a stationary random field A.) As usual, the homogenized coefficients are then given by 






6 


E 




They can be approximated by symmetric or non-symmetric approximate homogenized coefficients: 


t. /ffi) 

? ^hom 


f = lim E 

Ttoo 




6 


= lim E 

Ttoo 


'ffiA(p)(v4l + o] ■ (2.7) 


We denote by ^ := E • A(p)(V(^J|^ -h 0 


the non-symmetric approximate homogenized coef¬ 


ficients. When ^ is fixed, we simply write (^t for for (/>j, etc. 
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Difference operators. The aim of this article is to understand how depends on p for p close 
to 0. We shall first study the easier map p i—)• seen as a function of the approximate corrector 

Following physicists we introduce for all n € N a difference operator 5’^”^ acting generically on 
measurable functions of (ff, and in particular on approximate correctors as follows: for all C N, 

:= . 

This operator yields a natural measure of the sensitivity of the corrector (j)if with respect to the 
perturbation of the medium at inclusion This is to be compared to the vertical derivative used 
in [14] in the context of quantitative stochastic homogenization and to the randomized derivatives 
introduced by Chatterjee [7] in the context of Stein’s method (see also the Hoeffding decompositions 
in [20], where these randomized derivatives are used up to any order). For all finite T C N, we 
further introduce the higher-order difference operator 6^ = HnsF > more explicitly, this difference 
operator 6^ acting on approximate correctors 4)fp (for any if C N) is defined as follows: 

|i"l 

:= E = E (2.8) 

1=0 gcf gcF 

\G\=l 

with the convention ■= (t>ip ■ Physicists have introduced such operators to derive 

cluster expansions (see ]32]), which are used as formal proxies for Taylor expansions with respect to 
the Bernoulli perturbation: up to order k in the parameter p, the cluster expansion for the perturbed 
corrector reads, for small p > 0, 

ni ,712 GE^p'^ 

distinct distinct 

which we rewrite in the more compact form 

E = E (2.9) 

j=0 pcEip) j=0 fce(p) GcF 

\F\^j \F\=j 

where X]|G|=j denotes the sum over j-uplets of integers (when j = 0, this sum reduces to the single 

term G = 0). Intuitively, is expected to be close to a series where terms of order j involve a 
correction due to the interaction of j inclusions (and therefore derivatives of order j). Whereas the 
cluster formula for in Corollary 2.2 below holds under the mildest statistical assumptions on 

the coefficients, the validity of the expansion (2.9) is expected to require strong mixing assumptions, 
cf. [25] for the first order. This illustrates again the fact that averaged quantities (e.g. homogenized 
coefficients) are better behaved than pathwise quantities (e.g. correctors). 

For convenience, we also set 

5^(j)^ := +i - x), (2.10) 

that is, in terms of gradients, 

1^1 

V5f</>|( = ^^(-1)1^1-' + 0 + (2.11) 

1=0 gcf Gcf 

|G|=i 

By the binomial formula Yl\=o (^T^) = 0i have: V5^4)if = V5^4)if for all F / 0 and 

ff C N, = V())^ -|- ^ for all 77 C N, and, for all finite sets F,G, H C N, 

= Y, V<5f 

SCF 


( 2 . 12 ) 
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Inclusion-exclusion formula. When the inclusions are disjoint, we have 

cijp) = c-W. (2.13) 

n£E(p) 

However, when inclusions may overlap, this formula no longer holds since intersections may be ac¬ 
counted for several times. In the rest of this subsection we define a suitable system of notation to 
deal with these intersections. 

For any (possibly infinite) subset Fi C N, we set Ae := A + Ce, where Ce ■= {A' — A)\j^ 
and Je '■= Note that J|„} = = J„, and For non-necessarily disjoint 

inclusions, is then given by the following general inclusion-exclusion formula: 

(j(v) — ^{n} ~ ^{ni,n2} + 112,113} — ■ ■ ■ 

ni<n2^E^'P'^ ni<n2<n:i^E^‘P^ 

00 

= E(-1)‘*' E Cf. (2.14) 

k=l fce(p) 

\F\=k 

Since the inclusions J^s have a diameter bounded by 2R and p{B 2 r) is almost surely finite, the 
sum (2.14) is locally finite almost surely. Recalling that by assumption (2.1) the degree of the 
intersections of the inclusions is bounded by F, we deduce that we must have Cj? = 0 for all \F\ > F. 
Therefore, the inclusion-exclusion formula (2.14) actually reads 

r 

C(P) = Y{-1)^^^ Ce. (2.15) 

k=l fce^p) 

\F\ = k 

We shall need further notation in the proofs. For all E,F gN, E ^ 0, we set Je\\e ■— (flnsE Jn) \ 
(UneF Jn) and Jfp := (U^ge Jn) \ (UnsF Jn), and then 

C'eIIF := (^'- *^|^ ■= ("4'- Jf ■ 

M II ||i< 

In particular, we have Ce\\iz = Ce, C^^ = , and = 0. For simplicity of notation (except in 

the proof of Lemma 4.1), we also set 6 * 011 ^ = 0 = C 0 . The inclusion-exclusion formula then yields 
for all G,iF C N, G / 0, 

= ^(-l)l^l+iCs, 

SCH 

Co = E 

ScH 

Cg\\h = J;(-1)'^'C'5ug. 

ScH 

We shall also use the symbols ~, > and < for =, >, < up to constants that only depend on R, 
F, d, and A. Subscripts are used to indicate additional dependence of the constants, e. g. means 
that the multiplicative constant depends on rj, next to R, F, d, and A. Throughout, we will denote 
by C any positive constant with C ~ 1, whose value may vary from line to line. 

2.3. Statements. Our main result asserts the analyticity of the map p 1 —)• corresponding to 

the perturbed coefficients. 

Theorem 2.1 (Analyticity of the homogenized coefficients). Under the assumptions of Subsection 2.1, 
the map p 1 —)• is analytic on [0,1] and there exists a constant 0 < c < 1 such that, for allpo G [0,1] 


(2.16) 

(2.17) 

(2.18) 
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and all —pQ A c < p < (1 — po) A c, 


OO j 

a{po+p) _ a{po) , -|q^ 

^hom ~ ^hom Z_^ „-| ^hom J 
3 = 1 

where the series converges, and where, for any j > 1, denotes the (well-defined) j-th derivative 

of the map p i—)• at pQ. □ 

Since our proof is constructive, we obtain formulas for the derivatives. These formulas involve two 
approximation arguments: the addition of a massive term T~^ in the corrector equation to deal with 
integrability issues at large distances, and a hardcore approximation of the point process to deal with 
integrability issues at short distances. 


Corollary 2.2 (Formulas for derivatives). Let the assumptions of Subsection 2.1 prevail. We can 
construct a sequence {po )0 of hardcore approximations of the stationary point process p in the following 
sense: for any 9 > 0 , pg is an ergodic stationary point process on such that pQ C p, pe{Q) < 9 
a.s., and po f p locally almost surely as 9 f oo. For any F,G C N, denote by Aq , {Cg)p\\Q, {Cg)F 
the coefficients A^,Cp\\q,Cf corresponding to pg in place of p, and further denote by 4 >Fg^ the 
approximate corrector fip ^ associated with the coefficients corresponding to pg in place of p. 

Then, for all k >1 and all po S [0,1], the k-th derivative A^°^^ at po satisfies the following three 
equivalent formulas, for all 




(po),fc 

horn 




k\ lim lim 

Ttoo 0too 


\F\=k Gcf 


■GiE(po)\F 


fc E E. {C„)p\g||g(V0S™ + «) 

Too TOO 


T,e,^ 




\F\=k '-GcF 


( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 


where the limits exist and where the sums are absolutely convergent for any fixed T,9 < oo (recall that 
^\F\=k stands for the sum running over all the k-uplets of distinct positive integers). 

Moreover, in the case when the point process p satisfies E[p(Q)®] < oo for all s > 1, then the limits 
in 9 as well as all subscripts 9 can be omitted in the above formulas (2.20)-(2.22). Finally, in the 
case k = 1, and under the additional assumption that p{Q) < 9 q a.s. for some fixed 9o > 0, we can 
pass to the limit in T inside the sum in (2.21); for all po G [0,1], 


? ■ = E+ 0. ci”> (V0f + 0 

n 


(2.23) 


where the sum is still absolutely convergent. 


□ 


Formula (2.22) is the rigorous version of the so-called cluster expansion formula formally used by 
physicists (see [32]) as well as in [4]: it compares the homogenized coefficients corresponding to the 
coefficients obtained with a finite number of perturbed inclusions, fn particular, the fc-th derivative 
of with respect to p is obtained by considering k perturbed inclusions. Note that these cluster 

expansion formulas can be rewritten as in (2.9) using the difference operators defined in (2.8): for all 
k > 1, 


^ = k\ lim lim y" E 

^ ^ Ttooetoo ^ 
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rp P'iP0^\ P' _ P'iPo')\ P' 

where 6^ now acts on the random variable ^ ■ Ag ' {0){'V4>rpg (0) +^). For A; = 1, it essentially 
coincides with the formula obtained by Mourrat in [25]: 


{■ = li- S *: [«■ + 0 - «■ A"'”’''"’ (v4«+ «>; 

n 

Also note that, in the particular case when the inclusions JnS are disjoint, formula (2.21) takes 
the following simpler form: for all po G [0,1] and A: > 1, 

°°|F|=fcneF 

which further reduces, under the additional assumption that E[p(Q)^] < oo for all s > 1, to 

= A:!lim V V E + 0 • 

Ttoo ^ ^ ^ ^ L s 5S 

\F\=kn^F 


As a direct consequence of Theorem 2.1 we obtain the following universality principle, well-known 
by physicists: at hrst order in the volume fraction of the perturbation, the perturbed homogenized 
coefficient does not depend on the underlying point process p. More precisely. 


Corollary 2.3 (First-order universality principle). On top of the assumptions of Subsection 2.1, 
assume that E[p(Q)^] < oo. Then, we may define the volume fraction of the perturbation by the limit 


Vn := lim 
^ Ltoo 


E [ILQ n UneE(p) '^"1] 

F 


(2.24) 


and there exists some matrix K such that for all p >0, 

■4ll=Al + K«'p + 0(4)- 

If the point process p is independent of A, of A' (or else o/(A(,^)„ in the particular example (2.2)) and 
of the random volumes \J°\’s, then the constant K does not depend on the choice of the underlying 
point process p. □ 

Since the formulas given by Corollary 2.2 for the A;-th derivative of A^^^ at 0 involve terms 

of the form E[^^^ fiQm, ■ ■ ■, Qnk)]^ they depend on moments of p up to order k, so that stronger 
dependence on the point process p is expected for higher-order terms (see indeed [32, p. 493-494]). 

Formula (2.23) for the first derivative has the advantage of being exact (there is no limit left in 
T), and, at po = 0, it is given by the solution of the corrector equation corresponding to a single 
inclusion. In particular, this makes explicit calculations possible for spherical inclusions, and allows 
us to prove the celebrated Clausius-Mossotti formula in a very general context. 


Corollary 2.4 (Electric Clausius-Mossotti formula). On top of the assumptions of Subsection 2.1, 
assume that the inclusions are spherical, i.e. Jn = Bji{qn), and that both the unperturbed and per¬ 
turbed coefficients are constant and isotropic: A = aid and A' = /3Id. Denoting by Vp the volume 
fraction (2.24) of the perturbation, we then have, for all p >0, 


A 


ip) _ 

horn 


a Id -|- 


ad{/3 — a) 
(3 + a{d — I) 


Id Up 0{vl). 


□ 


As pointed out in the introduction, all our results also hold for linear elasticity. This allows ns to 
give the first rigorous proof of the elastic Clausius-Mossotti formula for random inclusions. Recall 
that an isotropic stiffness tensor A has the form : A : ^ = -|- ^(Tr(^)^, where G and A are the 

Lame coefficients, to which we associate the bulk modulus K = X-\- 2Gjd and shear modulus G. 
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Corollary 2.5 (Elastic Clausius-Mossotti formula). On top of the assumptions of Subsection 2.1, 
assume that the inclusions are spherical, i.e. = Bji{qn), and that both the unperturbed and perturbed 
stiffness tensors A and A' are constant and isotropic, and denote by K,G > 0 and K',G' > 0 their 
respective bulk and shear moduli. Let be the stiffness matrix of an isotropic medium of bulk 
modulus Ki = K + {K' — and shear modulus Gi = G + {G' — G) > where we have set 


d'^K + 2{d + l)id-2)G 
" “ ^ 2d{K + 2G) 




2G 


d-1 

d 


(2.25) 


Denoting by Vp the volume fraction (2.24) of the perturbation, we then have, for all p >0, 

= A + A^^Kp + 0{vl). 


□ 


Corollaries 2.4 and 2.5 treat spherical inclusions, in which case the solution of the corrector equation 
with a single inclusion can be calculated explicitly, so that (2.23) can be turned into an explicit 
formula. In the case of ellipsoidal inclusions, explicit calculations can also be made in terms of the 
so-called depolarization coefficients in the electric case (see [31]), or in terms of the Eshelby tensor 
in the elastic case (see [ 10 ], and also [26] for more precise analytic computations), so that an explicit 
formula for the first derivative can also be derived. The comparison of these results for spherical 

and ellipsoidal inclusions illustrates the fact that the first derivative already heavily depends on the 
geometry of the microstructure (see e.g. [32, Section 19.1.2]). 

An explicit formula could in principle also be obtained for the second derivative at po = 0 for 
spherical inclusions, since the corrector equation for two disjoint spheres can be solved analytically 
as well (see [30] and [18, Section 5]). 

Formulas (2.20), (2.21) and (2.22) for the derivatives as given by Theorem 2.1 are expressed as 
limits in terms of the approximate corrector gradient. For practical purposes, it may be important 
to prove rates of convergence for these limits. This is a quantitative ergodic result and therefore 
requires quantitative assumptions. In what follows we assume that a quantitative convergence result 
is available for the convergence of := E[^ • + ^)] to (through the convergence of 

to V(/>^^^) and show how this rate is inherited by their derivatives with respect to p. 


Corollary 2.6. On top of the assumptions of Subsection 2.1, assume that E[/ 9 (Q)^] < oo for all s > 1, 
and further assume that there exists a function 7 such that, for all T > 0 and p E [0,1], 


E[lV((/)P-4(J)|2]<7(r)2. 


(2.26) 


Let p G [0,1] he fixed. Recall the formulas for the approximate derivatives of A^^^: for all k > 0, 


\F\=kGcF 

Then, there is a constant C ~ 1 such that, for all A: > 0, we have 


(2.27) 


i(p)A _ ^(p)A 


2T 


< klG’^-tiTf 


In particular, if ^{T) for some a > 0, this yields for some constant G 1; 


A. 


(p)A 


-A 


(p)A 

horn 


< klC'^T- 


□ 


It is not clear to us whether Corollary 2.6 is optimal, and symmetric approximations could yield 
better rates. Such improvements, which would require nontrivial arguments based on quantitative 
homogenization theory, are not the goal of this article. 












12 


M. DUERINCKX AND A. GLORIA 


The optimal expected rate 7 for the approximate corrector gradient in (2.26) is as follows: 


7(r)2 


'T-i, ifd = 2; 

T-3/2, ifd = 3; 

T-2logT, ifd = 4; 

_r-2, ifd>4. 


It holds for instance under the assumption that the coefficients satisfy a spectral gap estimate (see 
[16]), and therefore covers the example of the Poisson point process {qn)n with inclusions J„ defined 
as the intersection of the Voronoi cells with the ball of radius 1 centered at Qn, with hxed conductivity 
A' = Ai in the inclusions, and T = ^2 in the reference medium. In higher dimensions, extrapolations 
techniques wrt T are needed (cf. [13]). 

In the following remark, we shortly discuss the need for assumption (2.1), that is hniteness of the 
degree of intersections between the inclusions. 


Remark 2.7. Assumption (2.1) is crucially used in the proof of the new family of improved energy 
estimates (Proposition 4.6), which are in turn at the core of the analyticity result. It is unclear to us 
how this assumption can be relaxed in general. In particular, this approach cannot treat the natural 
example of a Poisson point process with spherical inclusions of hxed radius (that is, p = {qn)n is 
a Poisson point process and Jn := B{qn), with (say) hxed conductivity A' = Ai in the inclusions 
and A = A 2 in the reference medium). In this specihc example, the random helds A^^i’s all satisfy 
a spectral gap estimate (in the form used in [12, 16]) and the point process satishes ]E[/3(Q)^] < 00 
for all s > 1 (we even have < 00 for all c > 0). The quantitative results of [12, 16] then 

provide additional analytical tools, as used in the discrete setting by Mourrat [25], which allow to 
prove a weaker version of Proposition 4.6 and conclude that the map p 1—)• is (at least) (7°° 

on [0,1] with derivatives given by the same analytical formulas as before. Yet, this is not enough 
to prove analyticity. In addition this approach makes quantitative assumptions on the random helds 
themselves (and not only on the point process), which contrasts dramatically with the results of this 
paper, and yields limitations in terms of applications (it is however enough to have the Clausius- 
Mossotti formulas in the form of Corollaries 2.4 and 2.5). Likewise, a quantitative approach (this 
time based on recent contributions by Lamacz, Neukamm and Otto [21]) allows to prove that, in the 
case of Bernoulli bond percolation on the integer lattice, the homogenized conductivity is (at least) 
C°° below the percolation threshold. We refer the readers to [9] for details on the results for the 
Poisson point process with sphrerical inclusions and for diffusion on the percolation cluster. □ 


Before we turn to the proofs, let us emphasize an observation by Anantharaman and Le Bris in 
[4] on the regularity of the corrector with respect to p in the case of disjoint inclusions — which 
dramatically contrasts with the analyticity of 

Remark 2.8. By testing the equation —V • — 4>) = V ■ + ^) in probability, we 

have 

- ^)\^] < E[]C'(^’)]2(1 + |V(^|2)] ~ p. 


We believe this scaling is optimal, so that the map [0,1] —>• L^(n) : p i—)• V(/>*'^^(0) is expected 
to be nowhere differentiable. Since we prove that the homogenized coefficients are analytic, this 
illustrates that averaged quantities behave much better than pointwise quantities (like the corrector 
gradient). □ 


3. Strategy of the proof 

In this section we present the strategy of the proof of Theorem 2.1. The key ingredient is a 
new family of energy estimates, the proof of which essentially combines combinatorial and induction 
arguments. When inclusions are disjoint, the combinatorics is signihcantly less involved than in the 
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general case of non-necessarily disjoint inclusions. In order to focus only on the core of the proof 
of Theorem 2.1 and to avoid additional combinatorial technicalities in this presentation, we shall 
momentarily assume that the inclusions are disjoint. 

Fix some direction ^ S |^| = 1. The aim of this paper is to investigate the difference 

{• (•411 - •4ho„)«, 

and express it as a convergent power series in the variable p around 0. Since the approximate correctors 
behave much better than the correctors themselves, we start with the analysis of the approximate 
difference 

:= f . (4> - Ar)(, 

for fixed T > 0. Indeed, the approximate differenee is a good proxy for the differenee since limy Ay — 
by (2.7). Next we rewrite the approximate difference in a form which is more suitable for the 
analysis. By definition, 

= E[e • + 0] - E[^ • AiV(j)T + 0] 

= E[e • + 0] + E[e • - ct>T)]. (3.1) 

The first term is already in a nice form (since it is of order p), while the second term is not (recall 
Remark 2.8: an energy estimate would only imply that it is of order y/p). In the following lemma, 
we make use of the corrector equation to unravel some cancellations. 

(p) 

Lemma 3.1. The approximate difference A)^ satisfies 

A^P = E[(V(/.t + 0 • C'(P) (V4^) + 0] ■ (3.2) 

Proof. Using that A^'i = A C^p\ the second term of (3.1) turns into 
E[e • AV{4>f - fir)] 

= + 0 • - fir)] - E[V<))^) • - fir)] 

= + 0 • A(p)v(4 ^) - fir)] - + 0 • - fir)] 

- ■ A^P\Vfi^p + 0] + E[V((.J’^ • A{V(pT + 0] + + C)]- 

By symmetry of the coefficients A and C^P\ reorganizing the terms yields 

E[e • - fir)] = -E[VfiT ■ A^P\V(l)f + 0] + E[V<))^^ • x4(V</.t + 0 ] 

+ E[V0t • + 0] • (3.3) 

The sum of the first two terms of the right-hand side of (3.3) coincides with the sum of the weak 
formulations in probability of the equations 

_ V . A^P^ +0 = 0 and - V • AiVfir + 0 = 0, 

tested with (pT and respectively, so that (3.3) reduces to 

E[C . AV{fiP - fir)] = E[VfiT • + 0], 

and (3.2) then follows from (3.1). □ 

Assuming that the inclusions are disjoint, we may use the inclusion-exclusion formula (2.15) in the 
elementary form of (2.13), so that (3.2) turns into 

A^ = ^E [{VfiT + 0 • + Olneuw] > 
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or alternatively, using the constraint n G to replace by (j)^ 

[(V</.T + 0 • + e)ln6E(P)’ • 

n 

Note that this sum is absolutely convergent since the C'^^'^’s are assumed to have disjoint supports. 
As l„gg(p) only depends on bn^ and {V(f>T + 0 ' {V(f)^ does not depend on bn \ we 

have by independence (2.3), using that bn'^ is a Bernoulli random variable of parameter p, 

= p ^ E [(V</.T + 0 • + o' • 

n 

We further decompose the right-hand side 

A^P =P^E [(V0r + 0 • + 0 

n 

+ p^E [(V 0 r + 0 • - 4 ">)] , 


and observe that the second sum is a difference of the same nature as , which begs for an induction 
argument, and the following lemma is indeed proved by induction (see Lemma 5.1 for a more general 
statement). 


Lemma 3.2. Assume that the inclusions Jn’s are disjoint and that E[/9(Q)®] < oo for all s > 1. For 
all k >0, all T > 0, and all p G [0,1], we have 

k 

fp) ^ Ai _L ^k+l T^{p),k+1 


i=i 

where, for all 0 < j < k, the approximate derivatives Alp and the error are given by 

-EE® • c^"^(v4+o], 

\F\=jneF 


Eip),k+1 ^ ^ E + 0 


|F|=fc+lneF 

and the sums in (3.5) and (3.6) are absolutely convergent. 


(3.4) 

(3.5) 

(3.6) 
□ 


Since the combinatorics in the proof of Lemma 5.1 is not more involved than for the proof of 
Lemma 3.2, we refer the reader to the proof of the former. 

If we can prove that \e!^'^’^\ < for all A: > 1 and for some constant C ~ 1 (independent of T > 0 
and of p G [0,1]), then we can easily pass to the limit T f oo in the expansion (3.4) and obtain a 
convergent power-series expansion for the exact difference A^^) itself around p = 0. 

The following lemma shows that a new family of energy estimates is needed to control the error 
terms. We display the proof of this lemma, which is significantly simpler than the corresponding 
proof in the general case of non-necessarily disjoint inclusions (see Proposition 5.2). 


Lemma 3.3. Assume that the inelusions J^’s are disjoint and that E[p((5)^] < oo for all s > 1. 
Then, there is a constant C ~ 1 (independent ofT, of p and of the moments of p) such that, for all 
k > 0, T > 0, and p G [0,1], the error defined in Lemma 5.1 satisfies 


|£;W,fc+l 


sE*: E 

1=0 L|G|=i 


\F\=k-j 

FnG=0 


k+1 


1=0 |G|=j 


(3.7) 


□ 
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Proof. Let k >0. First rewrite the error as follows: 


^{p),k+i = ^ ^ E <fT ■ + f) 


\F\=kn^F 

Recalling identity + C for H CN, we deduce 




EMuFUin} 


+c = 


Urp , 


GCF 


GCF 


SO that (3.8) turns into 


Eip)^k+i = ^ ^ ^ E [v<5f f)T ■ 


|F|=fcGcFn^F 


or equivalently 


p(p),fe+i 


^ ^ ^ ^ E ■ Cfo>(V(5f 

j=0\G\=jn^G \F\=k-j 

Fn(GU{n})=0 


(3.8) 


(3.9) 


For all n ^ G and all maps /, we obviously have (compare with the more general statement (4.22)) 
^ f{F,G,n)= ^ f{F,G,n)- Y1 fiFU{n},G,n), 


\F\=k-j 
Fn(GU{n}) = 0 


\F\=k-j 

FnG^0 


\F\=k-j-l 
Fn(GU{n}) = 0 


so that we may rearrange the terms in (3.9) as follows: 




i=0|G|=jn^G ■- 

+ EEEe 

j=0 \G\=jn^G 


\F\=k-j 

FnG =0 


E 


FUGUfn} 


(fr 


\F\^k-j-l 
Fn(GU{n}) = 0 




By Young’s inequality and the fact that the inclusions J„’s are disjoint, this yields 
k 

^E E E V,5f^^0T|1 +IE[|V5f4^V; 

,=0 |G|=, 


+ EEE4 

j=0\G\=jnfG ^ 
and the announced result already follows. 


E 


FUGUfn} 


4>T 


\F\=k-j-l 
Fn(GU{n}) = 0 




□ 


In view of (3.7) it is enough to prove the following family of energy estimates: there exists G ~ 1 
such that for edl k > j >0 we have 


E 


E I E Vi 

|G|=i \F\^k-j 

FnG =0 


fugAp) 


< G^. 


(3.10) 


On the one hand, a straightforward energy estimate directly yields (cf. Lemma 4.5) 


E 


Evi 


{n} Ap) 

i YT 


< 1 . 


(3.11) 
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On the other hand, an induction argument yields for some C ~ 1 and all j > 0 (cf. Lemma 4.4) 


E 


\F\=3 


F 4 p )\2 


< CL 


(3.12) 


For j <2, this estimate already appears in [4] (with however the massive term approximation replaced 
by the approximation by periodization). As mentioned in the introduction, in view of Lemmas 3.2 
and 3.3, these uniform bounds (combined with the fact that the estimates are independent of p and 
combined with some invariance argument due to the structure of Bernoulli random variables, see 
Step 3 of the proof of Theorem 2.1 in Section 5.2) imply that p i—)• is C^’^ on [0,1]. 

Before we describe the complete induction strategy used in Section 4.3 to prove (3.10), let us start 
by showing it in action, proving the result for k = 2 based on the corresponding result for k = 1 
(that is, (3.11) and (3.12) for j = 1). This proof is instructive in three respects: it implements the 
general induction strategy in the first nontrivial step, it shows that we need to use several forms of 
the equation satisfied by (pT, and it suggests that the proof of these equivalent forms relies on 

combinatorial arguments. 


Lemma 3.4. Assume that the inclusions J^’s are disjoint and that E[/9((5)^] < oo for all s > 1. 
Then, for all T > 0 and p G [0,1], 


E 






m^n 


E 


^1 

n m^m^n 

E I Y 

n^m 


< 1 , 


< 1 , 

< 1 . 


(3.13) 

(3.14) 

(3.15) 

□ 


Proof. For notational convenience, we consider p = 0 only. In what follows we take for granted 
that the series we consider are all absolutely converging, which is indeed ensured for fixed T by the 
(suboptimal) estimates of Lemma 4.2. We split the proof into four steps. In the first step we give 
three forms of the equation satisfied by (pT■ In the second step we prove (3.13) based on one 

equation and (3.11). In the third step we prove (3.14) based on another equation, (3.11), (3.12) for 
j = 1, and (3.13). In the last step we prove (3.15) based on a third form of the equation, (3.11), 
(3.12) for j = 1, (3.13), and (3.14). 

Step 1. Equations satisfied by 6^'^’'^^(pT- 

Let m 7 ^ re. The equation satisfied by the difference cpT can be written in several forms, with 

perturbed or unperturbed operators. With the unperturbed operator, we have 

- V • 

= V . + 0 - V • C^^^VcpY + 0 - V • + 0- 

By the inclusion-exclusion formula in the simple form due to disjointness of 

the inclusions, the equation takes the form 

- V • = V • ^ (3-16) 

This equation will be used to prove (3.15). A combinatorial argument (which is elementary here 
because the difference operators are of order 2 only and the inclusions are disjoint) allows one to turn 
the equation satisfied by (prp into 


(3.17) 
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which involves a perturbed operator (with the partially perturbed coefficients ), and will be used 
to prove (3.14). The third and last version of the equation takes the form: 

- V • = V • + V • (3.18) 

which involves the completely perturbed operator, and will be used to prove (3.13). 

Step 2 . Proof of (3.13). 

The starting point is (3.18), the right-hand side of which only involves hrst-order differences of the 
unperturbed corrector (pT- Although the argument of the expectation in the left-hand side of (3.13) 
is stationary, the equation (3.18) is not stationary. We shall first obtain energy estimates associated 
with (3.18) which are localized in space, ft is only after summing these estimates over n and m, taking 
the expectation, and passing to the limit in the localization parameter that the desired estimate (3.13) 
in expectation will come out in the form 


E 


n^m 


< CE 




(3.19) 


to be combined with (3.11). 

For all > 0, we then introduce a cut-off function XN for B]\f in B 2 N such that |Vx 7 v| ^ 1/A^i 
and test equation (3.18) with test function (pT S This yields for all n 7 ^ m after 

integration by parts, using the properties of XN, and rearranging the terms: 


J Bn B2N 

JB2N •tB2N 

We use Young’s inequality on each term (to ultimately absorb part of the right-hand side into the 
left-hand side), sum this inequality over n,m G N with n ^ m, and take the expectation to obtain 
(for a possibly larger C) 


[ E 


< c f E 


'Bn 

^ n^m ^ 

J B2N 

^ n^m 


+ 


— f E V + — f ^ 

^ Jb2n -I ^ JB2N 


where all the terms make sense and are hnite by Lemma 4.2. Since all the arguments of the expec¬ 
tations are now stationary, one may get rid of the integrals, which allows one to absorb the second 
right-hand side term into the left-hand side (choosing C > 0 big enough) and obtain 


E 


n^m 


< CE 


n^m 


c 


E 


n^m 


We are now in position to conclude: by taking f 00 we get rid of the second term of the right-hand 
side, so that (3.19) follows. 


Step 3. Proof of (3.14). 

The desired estimate is a consequence of (3.11), of (3.12) for j = 1, of (3.13), and of 


E 


^1 ^ 

*" n m,m^n 

< E 


n^m 


-fe 




(3.20) 
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The starting point is equation (3.18) that we first sum over m for m ^ n: 

1 v -Y 

m^m^n m,m^n 

= V- Y + V • Y 

m,m^n m,m^n 

Following the approach of Step 2, we test this equation in space with and the 

same cut-off xn- We obtain after summing the estimate over n, taking the expectation, and passing 
to the limit iV f oo, 


E 


^|v ^ 

*" n m,m^n 


< E 


El E 


m} 


n m,my^n 


-bE 


n m,m^n 


and hence, using that = 0 for n ^ mhy the disjointness of the inclusions, 

E' ■ ' ■ 


Ejv E 

n m,m^n 


<E +E 

^ n m^m^n ^ n m^m^n 

By the decomposition (pT + and the inequality < 1, the first 

right-hand side term turns into 


E 




m} 12 


n m,m^n 


< E 


E1^'*' 




-bE 


m^n 


The desired inequality (3.20) then follows from transforming the second right-hand side term as 
follows: we complete the sum over m, use the triangle inequality and the inequality 1 < 1 , so 

that 


E 


E".'"! E viW,#T 

m,m^n 


< E 


^ V(5^(/>r + E ^ I (/>r I 


(3.21) 


Step 4. Proof of (3.15). 

The desired estimate is a consequence of (3.11), of (3.12) for j = 1, of (3.14), and of 


E 


n^m 

< E 


-bE 


^ V(5^0T + E ^ II 


^1 ^ 

■ n m,m^n 

The starting point is equation (3.16), that we sum over n ^ m: 

1 Y - V • 2 IV ^ = 2V 'Y 

n^m n^m n^m 

Proceeding again as in Step 2, this yields 

2 


n} 


E 


n^m 


< E 


n^m 


(3.22) 


(Note that since each term of the equation is stationary after summation over n and m, this coincides 
with the energy estimate in probability.) Since the inclusions are disjoint, we are left with 


E 


n^m 


< E 


Y^Jr.\ E 

n m,m^n 


(3.23) 
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By the decomposition = 6^"^^ cj)T + cj)T, this turns into 


E 


n^m 


< E 


n m,m^n 


+ E 


n m,m^n 


and the desired inequality (3.22) follows from (3.21). 


□ 


This lemma implies that p i—)• is on [0,1]. As already mentioned, this lemma illustrates 

the induction argument we shall use in the proofs of Lemma 4.4 and of Proposition 4.6. 

In the proof of Lemma 4.4, we shall always consider the equation for d^(j)T with coefficients , so 
that the right-hand side will only involve unperturbed correctors, and then sum over F the resulting 
energy estimate (hrst localized in space). 

The proof of Proposition 4.6 is more involved. Call P{j,k) property (3.10). We make a first 
induction on k and then on j. Note that at step k there are k different forms of the equation satisfied 
by 5^(pT (for \F\ = k). By Lemma 4.4, P{k, k) holds for all A: G N. Then, given P{k -|- 1, A: -|- 1) and 
P{i, 1) for all z < Z < A:, we shall prove P{k + 1 — j,k + 1) iteratively starting with j = 1. Indeed, 
P{k + 1 — j,k + 1) will follow from P{k + 1 — j',k + 1) for j' < j and P{f ,1) for all j' < I < k, using 
the form of the equation where the coefficients are A;-|-1 — J times perturbed. The last step P{0, A:-|-1) 
is similar to Step 4 in the proof above and relies on the equation with the unperturbed operator. 

To be more precise, in the case of disjoint inclusions, the family of equations is as follows (see 
Lemma 4.1 for non-necessarily disjoint inclusions): for all disjoint subsets F,G,H C N, with F,G 
finite, F U G 0, 




n£F n£G 


4. Auxiliary results and improved energy estimates 

4.1. Perturbed corrector equations. We start by making precise the equations satisfied by the 
map 6 ^for disjoint subsets T, G C N, which will be used abundantly in this paper. The proof of 
this lemma (like many other auxiliary results of this paper) is purely combinatorial. 

Lemma 4.1. For all disjoint subsets F,F[ C N, with F finite, F 0, and for all T > 0, the map 
defined in (2.8) satisfies the following two equations (weakly) in 

SCF 

and 

^ ^(_l)|5|+lv.C5||HV(fP''<(.|^^. (4.2) 

ScF 

More generally, for all disjoint subsets F,G,H C N, with F, G finite, F U G 0, and for all T > 0, 
the map defined in (2.8) satisfies the following equation (weakly) in 

- V ■ 4,11. (4.3) 

SCFUCG 

Proof. Let T > 0 be fixed. Without loss of generality we may assume that H = 0. We first 
prove (4.2), from which we shall then deduce (4.3). Equation (4.1) is a particular case of (4.3) with 
G = 0. 

Step 1. Proof of (4.2). 

Let F C N be a finite nonempty subset. By definition (2.11) of and the inclusion-exclusion 
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identity (2.16), we have 

If </.T - V • AV6(c^t=Y1 V • + 0 

HCF 

HCFSCH 

SCF HCF\S 

Recognizing the definition of 4>^, this yields 

())T - V • AVf ^t=Y • CsVsPYt, 

ScF 

and proves the validity of equation (4.2). 

Step 2 . A combinatorial identity. 

For any finite subsets K,L,M C N (with K and L non-empty), we use the following notation: 

^L\\M ■— ~ f U f fl ^ f U ■ 

^neK ^ ^neL ^ ^neM ^ 

In this proof, and in this proof only, when K or L is empty, we further set and 

'^L||M ~ check the following general, purely combinatorial identity: for any finite disjoint 

subsets [/, F C N and for any SCF, 

HOF\S 

It is obviously enough to prove this identity for U = 5 = 0. Setting G := F \ S, we need to prove 
that, for any finite subset G C N, 

(-1)I^ICg= J](-l)l^\^lG|j^^^. (4.5) 

HcG 

Using the inclusion-exclusion identity (2.17) in form of have 

E(-i)'“'"'C|fj"= E(-i)'“'"' E (-i)'®'*‘csii„ = E(-i)''’'+‘ E (-i)'°''''csiiH. 

HcG HcG ScG\H ScG HcG\S 

Using then (2.18) in form of Cs\\h = ^ 5^0 this turns into 

E(-i)'°'"'cJ=j" = E(-i)'"'*' E (-1)'“'"'Et-iU'c^suc/ 

HCG SCG HCG\S UCH 

= ^(_l)|S|+l ^ (-1)\G\{HUU)\_ 

|CG UCG\S HcG\{SUU) 

Using twice the binomial identity in the form ~ ^K= 0 , this reduces to 

^(_l)|G\/I|c«GW= ^(_l)|S|+l ^ (-l)l^l + l^lG5uC/lf/=G\S 

HCG SgG UcG\S 

S^0 ' 

= (-i)I^Icg (-1)1^1+^ = (-i)I^Icg, 

SCG 

57^0 


and identity (4.5) is proven. 
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Step 3. Proof of (4.3). 

Let F, G C N be two fixed disjoint finite subsets, with F U G ^ 0. Equation (4.3) (with H = 0) is 
obviously a direct corollary of (4.2) (with F replaced by E" U G and with H = 0) provided we prove 
the identity 

-G^V5f + Af,g = Bf,g, (4.6) 

where we have defined 

AF.a-.= E 

SCFUG 

ScFUcG 

ScFUcG 


Let us first rewrite Af,g in n niore suitable way. We appeal to the definition (2.11) of 

then make the change of variables FI 0 S FI and U UW W, and conclude by using (2.16): 

ScFUcG HcF\SWcG\U 

HcFWcG ScHUcW 

^ E (-1 )I^\^I+I''\^Ig^^^(V<^|(^^ + 0 - (4.7) 

HcFWcG 

We now treat Bf,g- The change of variables F \ S S yields 


Bf,g= J;(-i)"'I+'J;(-i)I'^\'"g, 


{F\S)UU\\S 


VS: 


SUiG\U)^U 


(4.8) 


UCG 


SCF 




UcG 


SCF 


UcG 


=--Bfg 


_. o 2 
~'^F,G 


We treat both terms B^q and BpQ separately. The combinatorial identity (4.4) and the change of 
variables H 0 S H yield 


bL=E(-i)'‘''"'E E 


UcG 


SCF HCF\S 




SUIC\U)^U^ 


UcG 


HCF 


SCH 


By the identity (2.12) in the form YlscH^4 >t ~ V5^'^^, this turns into 




UcG 


HCF 


and thus, by definition (2.8)-(2.10) of 5^'^^ and the change of variables U 0W W, 

UCG HCF WCG\U 

VFcG UCVF 
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Noting that, by the usual inclusion-exclusion formula, 




ucw 


we conclude that 


^F,G — ~ 


ucw 

U^0 

- w - 




HCF 


VECG 


For the second term BpQ in (4.8), we argue as in Step 1, and obtain 


HcF ucG WcG\U 

U^0 ' 

HcF WCG ucw 

Uj^0 

= E ^{-1)^Uh\+\g\w\cW(^^^wuh 

HCFWCG 

Combining (4.7), (4.8), (4.9) and (4.10) then yields 

Af.o - Bf.c = E E - C | 15 ) + {). 

HCFWCG 


which proves (4.6) by definition (2.8)-(2.10) of d^'^^4>T- 


(4.9) 


(4.10) 


□ 


4.2. Basic energy estimates. The advantage of the massive term approximations (pp is to localize 
the dependence with respect to the coefficients to a ball of radius y/T (up to exponentially small 
corrections). While this regularization in T allows us to get rid of convergence issues at infinity, 
convergence problems may also occur at short distances because of high concentrations of the point 
process p. In order to avoid such issues, we further assume that p has all its moments finite. The 
assumption that T < oo and the finite moments assumption are crucial to make rigorous all subsequent 
formal computations. Under these assumptions we shall prove some basic energy estimates that are 
uniform with respect to the regularization parameter T and to the moments bounds on p; these 
estimates will be substantially improved in next section. 


Lemma 4.2. Assume that E[p(Q)®] < oo for all s > 1. Then, for all L ^ 1, k > 0, H G s > 1, 
and T > 0, the following estimate holds 



S-| 


< Cf^piBnY’^] < oo. 


for some constant Ct ~t 1; where is as in (2.10), and where we write {f)L{x) '■= 

for the local quadratic average of any map f. □ 


Proof. Since our argument is deterministic (we take the expectation only at the very end), we can 
assume w.l.o.g. that H = 0. By (4.1) in Lemma 4.1, (fx satisfies 

0T - V • ^t=Y1 • Cs\\F\s^5U^(t>T. 

SCF 


Let z G and set pf^{x) := e ^Jth c > 0 to be chosen later. Testing this equation with 

(j)T in the whole space, and noting that |Vp|,| = cpp/y/T, we obtain the starting point for a 
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Caccioppoli’s inequality 

— f + f Vt \'^^ ^ 

J JKd JR‘i tri 


_ Vt\^S!<Pt\\V6P^cI>t\ 

ScF“'-^s||f\s 

+ ^ E / ^t\S!M |v<5f / ^Tl^f </>r| |vf 

V j 17 J JfjnF\S y/T jRd 


SCF'' '^S\\F\S 

It is crucial to note here that, for fixed F, the sets Js\\f\Si S G F, are all disjoint. By Young’s 
inequality, and choosing c > 0 small enough so that one may absorb all the terms but two in the 
left-hand side, this turns into 


T jRd ^ JRd ^ Fri 


ScF-'-^s\\f\s 


For L ~ 1, taking the square root of both sides yields 
{T-ks!<pT)Liz) + iV5^cPT)L{z)) < ( 


ScF-'-^s\\f\s 


< 


E 


ScF VJs\\F\s 








(4.11) 


Now note that, renumbering the sum in terms of F’\ S', and using that Js\\f\s — ^ whenever S = 0 , 
we get 


E E (/ < E E (/ ^Tiv<5f0TN , 

|F|=fcScF V-''^S||F\s / |F|<fc-l|S|<fc V’^'^SIIF / 

and hence, as C Jn C Bji{qn) for any n G S, 

EEf/ 'fTlV<5f''’<ATp) ’ < E zf/ fTlVjf.i.Tl^y E (4.12) 

7T^l_7. C,—Z? \ '^’Jf!\\F\F! } irnl^;. i n \'^^R(Qn) } \Q\^u 


\F\=kS(lF V'^S\\F\S 


\F\<k-l n Y^R 


\S\<k 

nGS 


We bound the last sum as follows: for any fixed n, recalling that by assumption (2.1) the intersections 
of the inclusions J„’s are of degree at most F ~ 1 , 


E! ^■^s||F^0 — E] — 




< 2 ’"-^ < 1 . 


(4.13) 


n^S 


n^S 


As we have t]^{x) < e z\IFt^ cRI\fT ^ ^ Bfi{qn), we can then deduce from (4.12) and 

(4.13): 

|F|=fcScF V-''^S||F\S / |F|<fc-l n 

We then sum (4.11) over 1^1 = k, k > 1, and use the above estimate to get for any z G W^, 

SY) := T--2 ^ (5f0r)F(^) + ^ (Y^tUz) 

\F\=k \F\=k 

c{L+R)/VT ^ ^-c\q„-z\/VT 


<gC(F+KVVF^ ^g-c|g„-.|/VF ^ {V6^^T)R{qn). 

\F\<k-l 
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Combining this with (2.5), we conclude by induction that, for some (deterministic) constant Ct ~t Ij 


e 

j=l 


-c\qn^ -z\Is/T -c\qni_^ -q^il/VT 


n 

i=2 


=: 4 ( 2 ) 

It only remains to compute the sum I^{z). For that purpose, we compare sums to integrals 


■jR/Vr 


4 (z) < e' 
and hence, 

4{z) < 


E [ ■■■[ 


^-c\xi-z\/VT 


ni,...,nj 


RiQni) 


j 


i=2 


''-'^l"^-"l/^p(i?R(xi)) n dxi... dXj. 

i=2 ^ ^ 


Taking expectation of Ii^{zY, for some s > 1, and applying the triangle and the Holder inequalities, 
we obtain 

r ^ 

E[4(z)"]^/* < e^^^/^E[p{BRY4/^ / ... dxj, 

which finally gives, by an obvious change of variables, 

E[4(^)"] - < e^^^/^E[p{BRYY~^ J = C^T^^/^e^^^/^E[p{BRYYK 

and the announced result is then proved. □ 

Based on this deterministic estimate, we prove a lemma which will be crucial to give sense to formal 
calculations, and implies the absolute convergence of all the series we will be considering for fixed T, 
at least under the additional assumption that moments of p are finite. Note that the result obviously 
also holds for V(/>y replaced e.g. by . 

Lemma 4.3. Assume that E[p((5)^] < 00 for all s > 1. For all T > 0, and k > 1, we have 

E E [IC’ncI |V<5f</> t| (1 + |V<)>^ 1 )] < 00 , (4.14) 

|F|=fcGcF 

^e[|Cg||v5P''</>4(i + |v</.|:|) 

|F|=fcGcF 

^ E [\Cr\ (1 + |V</.t|)( 1 + |v</.^|)] < 00 . 

|F|=fc 


< 00 , 


(4.15) 

(4.16) 


Proof. We only prove (4.14); the proofs of the other statements are similar. Let /c > 0 be fixed. By 
stationarity we add a local average over the ball Br, say, with L ~ 1, we apply the Cauchy-Schwarz 
inequality, and note that, for all x G Br, 

\Ch{x)\ < la;eJn,VneLI < lxeSfl(q„),VneF = lg„eFH(a:),VneF < IgneBfl+i.VneF =: Xl{H), 

SO that we can write by the change of variables F T U G, 


4 < E E 

|F|=fcGcF 

s E E*^ 

|F|<fc |G|<fe 


XL{F\G){V5ff>T)L{l + {V(^^)L 
XF(T)(V<5f0r)F(l + (V0f:u«)i) 
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Using the deterministic estimate (2.5) in the form of ^ this yields 


< T2]E 


\F\<k ^ ^ |G|<fc 


The first sum can be estimated as follows: since (”) < n*/i! < {enjif and ^ t) 


E = E 


< N Y < 


‘-gn&Bii+L,VneF ^ \ i J 

|F|<fc |F|<fc i=0 ^ ^ 


E 


p{Br+l)^, 


so that we conclude by Lemma 4.2 and the Cauchy-Schwarz inequality that 


gk < r^/2]E[p(i?^^^)2fc]l/2E / ^ {V5^c^t)l 

-^\G\<k 


1/2 


< oo. 


□ 


We now turn to energy estimates that hold uniformly with respect to T and the moments bounds 
on p. The following two estimates will be further improved in the following section. 

Lemma 4.4. Assume that E[p(Q)^] < oo for all s > 1. Then, there exists a constant C ~ 1 
(independent of T and of the moments of p) such that, for all k >0 and T > 0, 


E 


\F\=k 




□ 


Proof. For A: = 0, the result E[|V(/)t + Cp] ^ 1 reduces to the basic energy estimate on the modified 
corrector (see (2.5)). We now argue by induction. Assume that the result holds for some fixed A: > 0. 
From (4.1) in Lemma 4.1, we learn that 5^4>t satisfies on 

cfT-V- <Pt=Y. • Cs\\F\s^5l^^(>T. 

SCF 


We test this equation with where XN is a cut-off function for Bjq in B 2 N such that |VxArl ^ 

\/N. This yields 


[ |Vf</.rP<E/ (4.17) 

J B]\f g(~p’'B2N 

+ 4 A |4f^,T| + 4 ^ |V«P0 t|. 

J Bom c’^t'JBom 


We then take the expectation, sum over \F\ = k — 1 (all the sums are convergent by Lemmas 4.2 
and 4.3), and divide by N^: 


E 


< E 

rv_/ 

/ E EUs.nJVif^TlIVifVl 


-^^\F\=k+l 


\F\=k+lScF J 








Since each sum above is absolutely convergent and dehnes an integrable stationary random field (the 
expectation of which obviously does not depend on the point it is taken), this inequality also takes 
















26 


M. DUERINCKX AND A. GLORIA 


the form 


E 


|E|=fc+l 


< E 


|F|=fc+l SCF 




^ \sl[<h-\\ysc<h\ 

|E|=fc+l 




E ED,„«iif0Tiivif'Vi 

|F|=fc+l SCF 


Taking the limit f oo then yields 


E 


|E|=fc+l 


< E 


|F|=fc+lScF 


(4.18) 


By Young’s inequality and the disjointness of the sets Js\\f\s^ S G F (for fixed F), in the form of 
Z^5CF l./s||F\s < 1 , we may absorb part of the right-hand side into the left-hand side, and obtain 


E 


< E 

rs_/ 



EE|=fc+l ^ 


L |F|=fc+l5cF 




< E 


5 


\F\<k |S|<fc+l 


where we used that J 0 ||j? = 0 . 

By assumption (2.1), proceeding as for (4.13), we have 

/p\ 

E D,„( 0 ) < ^ 1 *( 0 ) < E ( ,■ S 2 '- < 1 , 

|S|<fc-ri |S|<fc+i j=i 

so that (4.19) finally turns into 


E 


\F\=k+l 


< E 


\F\<k 


from which the desired conclusion follows by the induction assumption. 


(4.19) 


(4.20) 


□ 


For sums of size k = 1, the following result is easily proven as an energy estimate in 

the probability space; for general k it also holds but the proof relies on a subtle induction and 
combinatorial argument, which is presented in next section. 


Lemma 4.5. Assume that E[p((5)^] < oo for all s > 1. Then, for all T > 0, we have (uniformly in 
T and in the moments of p) 


E 


Ev<5E</>t 


< 1 . 


□ 


Proof. By Lemma 4.2 for k = 1, the sum is well-defined in F[(^^{W^) and satishes the 

following equation on 


f E - V • E = V • e +0- 

n n n 


We then test this equation with XNi'^n f’T) for some cut-off XN foi' Bj\i in B 2 N such that |Vxa| ^ 
l/N. Since is stationary, we may proceed as for the proof of (4.18) in Lemma 4.4, and 
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obtain after taking the expectation and the limit N oo (or equivalently testing the equation in 
probability): 


E 




< E 




< 

rs_/ 


E 


n ^ n ' 


By assumption (2.1), proceeding as for (4.13), we have X^„lj„(0) < 1, so that, using in addition the 
decomposition 1 + ^ (1 + |V(/>tP) + we obtain 


E 






{n} 


<pT 


< 


E[l + |V((.Tp]+E 




< 1 , 


where the last inequality follows from Lemma 4.4 with k = 1. 


□ 


4.3. Improved energy estimates. In this section, we prove the following generalization of Lemma 4.5 
to any order k in the following form (choosing j = k in (4.21) below): there is a constant C ~ 1 such 
that, for all A: > 1 and T > 0, 


E 


V6^(j)T 

\F\=k 




and we give an interpolation result between this inequality and the energy estimates of Lemma 4.4, 
that is, a , L^(n))-type estimate. 


Proposition 4.6. Assume that E[/9((5)*] < oo for all s > 1. Then, there exists a constant C ~ 1 
(independent ofT and of the moments of p) such that, for all T > 0, k > 0, and 0 < j < k, 


Sf := E 


E 


E 




\G\=k-j \F\^i 

Fna=0 


< c 


k+1 


(4.21) 


□ 


Proof. We proceed by induction and split the proof into two steps. 


Step 1. Preliminary. 

For k = j = 0, the estimate Sq = E[|V(/)'r + < 1 reduces to the energy estimate for the modified 

corrector. (Note also that the estimate 5^ < 1 already follows from Lemma 4.5 — but this will not 
be used here.) We now argue by a double induction argument. Since the result is proven for k = 0, 
we may indeed argue by induction on k: we assume that Sj < for all 0 < j < fc' and for all 

0 < k' < k, and shall prove that for all 0 < j < A: + 1. Since Lemma 4.4 implies the 

desired result for j = 0, we may as well argue by induction on j: we further assume that < 0^+“^ 

for all 0 < j' < j, for some 0 < j < A: + 1, and shall prove that <S'j+i < 

Before we turn to Step 2, we state another combinatorial inequality we shall need in the proof. Let 
G, S' C N be finite fixed disjoint subsets. We claim that 


We first rewrite 


E 

\F\=k 

Fn(GuS)=0 




(fr 


< 


| 5 | 

EE 

/=0 \L\=l 
LCS 


\F\^k-l 

Fn(GUL )=0 


^G,S 


E 




4>t = Y/ lFnS=0V(5|^‘^(/)T 


FUG, 


\F\=k 

Fn(GuS)=0 


\F\^k 

FnG=0 


where we can decompose, by the usual inclusion-exclusion argument. 


|5| |5| 

lFnS=0 = 1 - ^Fns^0 = 1 - ^(-1)*"^^ Y^ Ifcf = Y.G 

l=l \L\=l 1=0 

LCS 


1)^ Y ^lcf, 

\L\=l 

LCS 


(4.22) 
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SO that Sq g becomes, by a change of variables, 

|S| |S| 

^G,5 = E(-1)'E E = E(-1)' E E 

1=0 \L\=l \F\=k Z=0 \L\=l \F\=k-l 

LCS FnG=0 LCS Fn(GUL)=0 


and the claim (4.22) then follows from the triangle inequality. 

Step 2. Bound on 

Let G C N be a finite subset. By Lemma 4.2 we may sum equation (4.3) of Lemma 4.1 (with 
H = 0) for over \F\ =j + l, F D G = 0, which yields the following equation for the sum 

X) |F|=i+i on 

FnG=0 

\F\=j + l \F\=j + l 

FnG=0 FnG=0 

\F\=j+i ScFUcG 

FnG =0 

UcG |S|<3+1 |r’|=j+i-|S| 

SnG=0 Fn(GuS )=0 

We test this equation with XAf X]|F|=j+i FnG=0 XN is a cut-off function for Bn in 
B 2 N such that |VxAr| < ^/N, we take the sum over |G| = {k + 1) — {j + 1) = k — j (which is again 
absolutely converging by Lemma 4.2), take the expectation, and then use stationarity to pass to the 
limit N 'I 00 , as in the proof of (4.18) in Lemma 4.4. With the notation a A 6 = min{a, 6}, this yields 


^ IE 

‘ E 

E 

E 

(-l)l"'+'‘'l+'CsuG||GW E 

2- 


I 

II 

UcG 

|5|<a+i)Ar 

SnG= 0 ,SuU ^0 

|F|-j+l-|5| 

Fn(GuS)=0 



where the additional restriction S 0 U jtz 0 follows from the fact that G^ut/ncy^/ vanishes identically 
otherwise and where we have further restricted to \S\ < L since by assumption (2.1) there is no 
intersection of degree larger than L. Since we have IG^uf/ycyf/l < Ijgl (using here notation 

\j^ = 1), and the J[/||G\c/’s are disjoint for [/ C G (for fixed G), we deduce 




E 


E^ 

\G\=k-jU(ZG 




G||G\G 


E 


Gs 


|S|<0+i)Ar 

SnG=0,SGUi^0 


E 




|F|=j+l-|S| 

Fn(GuS )=0 


As for (4.20), we have X^| 5 |<j+i Ijs(O) < 1, so that by the Cauchy-Schwarz inequality, this estimate 
turns into 


s’;xi < E 


E E E 


Gs 


\G\=k-j UcG 


|s|<(i+i)Ar 

SnG=0,SUG?^0 


E 




|F|=j + l-|S| 
Fn(GuS)=0 


2n 


Now using the decomposition V= XIrcS (that is, (2.11) with H = 0, 

G F U (G \ G) and F S), together with the observation that 


Gs 


E “^2) - (E ~ E 

RCS ^ ^ RcS ^ Res 
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which follows again from combining the Cauchy-Schwarz inequality with inequality 1 ^ 1, 

we obtain 




< E 

r\j 

W 

^Ju 

E > 

Js E 

E 

Vhfu^u(G\u)^^ 

2- 



L|G|=fc- 

jUcG 

|5|<0'+i)Ar 

SnG= 0 ,SUU ^0 

Res 

\F\=j + l-\S\ 

Fn(Gu5)=0 



(j+i)Ar 







2- 

< E E 

' E 1 

Ju ^ijk 

E > 

Js E 




j=0 

^\U\<k-j 

\G\<k-j 

Gnu =0 

|5|=i 

SnG =0 

ReS 

\F\=j+l-i 

FniGUUUS )==0 




where we have set <5^^ = 0 when simultaneously |G| = k — j and i = 0, and <5^^ = 1 otherwise. By 
(4.22) and the inequality 1 ^ 1, for any R G S and any G CiU = 0 = S Ci G, we have 


T^JsnJu 


2 

< 

/i+i-'i 

E 



\F\^j + l-i 

Fn{GUUUS)^0 


^ 1=0 

\L\=l 

LCUUS\R 

\F+j+l-i-l 

Fn(LURUG)=0 


j+l-i 

^ Y 1 

/=0 \ M =1 

LCUUS\R 


E 


^ ^FULURUG 


4>t 


Fn(LURUG)=0 


2 


(4.23) 


and hence we obtain, using 'F\u\<k^Ju ~ ^ again, and using U L) {S \R) GN\{GL) R), 


(j+l)AF j+l-i 


rifc+l < \ ^ 

^j+1 ~ 

E E 

E 

^Ju E 



i=0 

/=o 

'-\U\<k- 

j \G\<k-j 





Gnu=0 


2- 




E 

ijs E E 





|5|=i 

Sr\G=0 

Res \L \=1 

LGUUS\R 

\F\=j + l-i-l 

Fr]{LURUG)=0 



(j+l)AF j+l-i 





2- 

S E 

E IE 

E 

E 

^Jr E 



i=0 

1=0 

^\G\<k-j 

|H|<i 

RnG—0 

\L\=l 

Ln(GUR)=0 

\F\=j + l-i-l 

Fn(LURUG)=0 



Successively using T < 1 in the form of E|R|<i ^ E|R|=i and El^Jl ^ we 

obtain by the change of variables LL) R R, 


(j+l)Ar j+l-i 





2- 

qk + 1 ^ \ ^ 

^j+1 ~ 

E E 

Eo, 

E 

E 



i=0 

1=0 

L|R|=i 

\L\=l 

\G\<k-j 

\F\=j + l-i-l 

- 



LnR=0 

Gn(LUR)=0 

Fn(LURUG)=0 



(j+l)Ar k-j 




2- 

^ E E*^ 

E 

E 



i=0 1=0 

L|R|=i 

\G\=l 

RnG=0 

\F\=j+l-i 

Fr\(GUR)=0 



or equivalently, recalling the dehnition of the and of the Sj's, 



k-j-1 

j+1 k-j 




2- 

qk + 1 

^j + 1 

< c«+i+i 

~ Z^ ^j+i 

+ EEi^ 

E ^Jr 

E 




1=0 

i=l 1=0 

L|R|=i 

\G\=l 

RnG =0 

\F\=j+l-i 

Fn(GUR)=0 



(4.24) 
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Using again the fact that ~ bound 


w 

1_ 

E 


2- 

= E 

M 



2- 

AR\=i 

\G\=l 

RnG =0 

\F\=j + i-i 

Fn{auR)=0 



■ |G|— i-\-l 

RCG 

\R\=i 

|F|=j+l-i 

FnG=0 



< 

rv-/ 


E 


E 

|G|— i-\-l 


|F|=j + l-i 
Fr\G —0 


2n 


ol+j+l 


SO that (4.24) turns into 


qk + l < 
^j+i ^ 


k—j—1 

E sSb' 


j+l k-j 

+EEs;:u. 


1=0 


i=l 1=0 


k j k+1 

E'Shi + E E 


1=0 i=0 l=j+l 


As the right-hand side only involves the Sji’s with k' < k or with k' = k + 1, f < j, we conclude 
that E by the induction assumption. □ 


5. Proofs of the main results 

In this section, we prove the analyticity of the perturbed coefficients (Theorem 2.1) and the analyt¬ 
ical formulas for the derivatives (Corollary 2.2), from which we further deduce the Clausius-Mossotti 
formulas (Corollaries 2.4 and 2.5). 


5.1. Approximate derivatives at p = 0. In this subsection we devise analytical formulas for the 
derivatives of the map p i—)■ at p = 0 under the assumptions that T > 0 and E[p(Q)^] < oo for 

all s > 1. We shall show in particular that is C°° at p = 0. These results, which rely on the 
improved energy estimates of Proposition 5.2, constitute the core of the proof of Theorem 2.1. 


Fix some direction ^ |^| = 1. As in Section 3, we consider the exact and approximate 

differences 

^(P) ^ ^ 

and we recall that limT’follows from (2.7). By Lemma 3.1 the approximate difference 
satisfies 

= E[(V0t + 0 • C(^’)(V0P + 0]- (5.1) 

By assumption (2.1), we may now appeal to the inclusion-exclusion formula in the form of (2.15), 
so that (5.1) turns into 


At^) = 


E( 

i=i 


\F\=j 


E 


{V^T + 0 • + mPcEM 


where the sum is absolutely convergent by (4.16) in Lemma 4.3. Using that the event [F C is 
by definition independent of the rest of the summand, and that we have i.i.d. Bernoulli variables of 
parameter p, this identity takes the form 



r 

^(-i)^+Ve 

i=i 


\F\=j 


(5.2) 
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which can be further decomposed as 

r 

= ^(-ly+V'E 


i=i 


J](v</>T + 6-c'u(v</-^ + e) 

\F\=j 


+ E 

i=i 




\F\=j 


where the sums are still absolutely convergent by (4.16) in Lemma 4.3. The hrst term of the first sum 
(i.e. corresponding to the choice j = 1 ) is of order p and coincides with the argument of the limit 
in (2.20) for k = 1. The second sum can be rewritten as a sum of errors of order at least p^, which 
can then be combined with the corresponding (higher-order) terms in the hrst sum, and an induction 
argument hnally allows us to prove the following decomposition: 

Lemma 5.1. Assume that E[p(Q)®] < oo for all s > 1. For any k > 0 and any p £ [0,1], we have 


P’^T 

j=l j=k+l 

where, for all j > k > 0, the approximate derivatives A^j, and the errors are given by 

4 E E (-i)"'\'"+'iE [vsfcfT • c^\giig(v<^|: + 0], 

|F|=iGcF 


fc+r 


i E-(p) 




^ ^(_l)inG| + lE + 0 


\F\=j gcf 
' ' •' |G|<fc 


and the sums in (5.4) and (5.5) are absolutely convergent for fixed T. 


(5.3) 

(5.4) 

(5.5) 

□ 


Proof. We proceed by induction. For k = 0, (5.3) reduces to (5.2). Assume now that (5.3) holds true 
for some /c > 0 . 

First of all, we decompose as follows: 


where the error reads 


G. 


ip)M _ 


^{p),k+l,k _ ^fc+1 Q{p),k 


Y, E (-1)"'\'"+'IE [vhf <).T • CF\G\\G^{cfY''^^ - </>r) 


(5.6) 


|F|=fc+l gcf 
' ' |G|<fc 

k+l 




|F|=A:+1 J=1 


GCF 

IGNj 


since the summand for G = T in (5.4) vanishes (cf. C' 0 ||^ = 0). 

Given \F\ = k + 1, recall that (4.1) in Lemma 4.1 (for H = 0) asserts that solves 


k-\-l 

^6^ cfT-V- cfT = E(-l)'^' E ^ 

i=i 


^G||F\GVhf'',/.T, 


GCF 

\G\=j 


and also recall that since _ ^F _|_ Gy^, ~ 4>t solves 
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Successively testing these equations with — and 6^4>T respectively (as for the proof of (4.18) 

in Lemma 4.4, still noting that all the sums converge absolutely by Lemmas 4.2 and 4.3), we get 



~f ^ L ^^ - (t^r) 

\F\=k+l 

|F|=fc+l 


Y VS^cl)T ■ A^V{4>§ 

\F\=k+l 


eMuf 


E(p)uF 


+ 0 



Hence, using the inclusion-exclusion formula (2.15) as before (cf. (5.2)) and the independence, this 
yields 


r^{p),k 

rp 


^{-ly+y ^ 


E 


i=i 


|G|=j 


GnF=0 




and hence, renumbering the sums. 


V 


fc+i 


riip),k 

\_T p 


k+V+1 

E (-ir'-y E 


E E 


j=k+2 


|i"l=i 


GCF 

|G|=fe+l 




(5.7) 


By the induction assumption (5.3) at order k and the decomposition (5.6), we thus have 


fc+1 fc+r 

, (p) _ ^\3 , 


=Y.P^K+p''^^Gt + Y. ^(-l)inG|+l]E Vdfcj)T-Cp\G\\Gi^^T ^ + 0 

i=i 


j=k+2 \F\=j Ggf 

' ' |G|<fc 


Combined with (5.7), this yields 
k+i fc+r+i 

A^^t=Yp^^t+ Y p’ Y Y (-1)''^'’^^1E 
1=1 


j=k+2 |F|=7 Ggf 

' ' |G|<fc+l 


V5'^<Pt-Cf\g\\g{^ 4>Y'^^ + i) 

and therefore (5.3) at step k + 1. 


that is = Ytl P’^T + Efcfc +2 P’^i 


□ 


We now prove that the approximate derivatives are bounded uniformly in T and in the moments 
of p, as a consequence of the improved energy estimates of Proposition 4.6. 

Proposition 5.2. Assume that E[p((5)®] < oo for all s > 1. Then, there is a constant (7 ~ 1 
(independent ofT and of the moments of p) such that, for any k>l, the approximate k-th derivative 
defined in (5.4) satisfies 

|A^| < (5.8) 

Likewise, for any j > k>l and any p £ [0,1], the error defined in (5.5) satisfies 

□ 

Proof. The estimates of the errors are obtained using the same arguments as for the estimates 

of the approximate derivatives Aj,^s, and we only display the proof of the latter. Since + '^ = 
'Esce'^^I^t (cf. (2.12) with G = iL = 0) and Cf\g\\g = Cf\g + Ei/cG,; 7 ^ 0 (-l)'^'C't/u(nG) 
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any G C F (cf. (2.18)), and C^wg = 0) we may rewrite formula (5.4) as follows: 

\F\=kGgF SCF 

'-V-' 

= ■ A*: 

■^T,i 

+ E E E (j)T ■ C'c/U(F\G)V'5f</’r] • 

|F|=fcGCFScF UdG 

' -V-' 

= :A^2 

We treat each term separately. By the change of variables G ^ G U U, S S D U, and F 
F \J G U S U U (with F, G, S, U disjoint), we rewrite ^ as 

^t,i = EEE E E E E 

1=0 i=0 j=0 \F\=k-l-i |G|=i-i |S|=i \u\=j 

GnF =0 sn(FuG )=0 un(FuGuS )=0 


SO that, by the triangle inequality. 


fc—1 k—l I 


|A 


T,1 


< 


EEE E E E ^ 

^=0 2=0 j=0 \F\=k-l-i |5|=i \u\=j 

SnF =0 un(FuS)=0 


ljFus|V5f''</>T| 


^ Vdf^^cPT 

\G\=l-j 

Gn{FuSuu )=0 


Recall from (4.22) in the proof of Proposition 4.6 that 



1^1 I5I 



sEEE E 


\GM-j 

GniFuSUU )=0 

u=0 S = 0 \W\=u |if|=s 
WCF HCS 

\G\ = l—j — u — s 

Gn{WUHUU)=0 


Hence, by the change of variables F F \ W and S S \ H, and the notation 6f,s,U,W,h = 1 if 
F, S, U, W, H are disjoint, and Sf,s,U,W,h = 0 otherwise, this yields 


k—l k—l I k—l—i i 

^r.ii^EEEEE E E E E ^^fsf,w,h 

1=0 i=0 j=0 u=0 s=0\F\=k-l-i-u\S\=i-s\U\=j\W\=u\H\=s 


X E 





\G\ = l—j — u — s 


Gn{wuHuu )=0 


We rearrange the sums suitably, and use the notation Ij,^ = 1 (so that we have = Ij^lj^^) to 

obtain 


k—l k—l I k—l—i 


II 


< 


EEEEEE E 

1=0 i=0 j=0 u=0 s=0 \U\ = j 

Hnu=0 


E 


^Jh 


E 


^Jf 


|F|=Al —/ —2 —22 


E D.IVil 


SUHUU 


4>t\ 


|,S|—i —s 
Sn(Huu )=0 


< 

E 

^Jw 


)] 


\W\=u 


\G\^l-j-u — s 



wn(Huu)=0 Gn(wuHuu)=0 
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Recalling that XIlcN (^) ~ ^ (4.20) (as a consequence of assumption (2.1)), we deduce from a 

multiple use of the Cauchy-Schwarz inequality 

k—1 k—l I k—l—i i 


K.iisEEE E E E E E 

1=0 i=0 j=0 u=0 s=0\u\ = i 

Hnu=0 


Jh 




SiJHyjU 




|5'|=i—s 
Sr\{H\JU)=0 


E 


^Jw 


\W\=u 

wn{Huu)=0 


E 


^^^GUWUHUU 


4>t 


\G\=l—j—u — s 
Gr\(WUHUU)=0 


and hence, by the Jensen inequality, 

k—l k—l I k—l—i i r 

K,.i<eEE E E*^ E E E 


■SUHUU 


(Pt\‘ 


1=0 i=0 j=0 u=0 s=0 \U\ = j 

Hnu=0 


|5'|=i—s 
Sn(HUU)=0 


'e 

E 

E 

^Jw 


2- 

^\u\=j 

\H\ = s 



\G\ = l—j — u—s 

-* 


k—l k—l I k—l—i i 

+ EEEEE‘ 

/=0 i =0 j =0 u =0 s =0 

Mnu=0 wn{Huu)=0 (jn{wu±iuu )=0 

By the changes of variables SUHUU'^U in the first term and WUHUU U in the second term, 
and using that E|iT|<fc Yl\w\<k ^Jh^Jw < this finally yields 


k-l j 


(5.10) 


E +EEe E E 

j=0 ^\u\=j J j=0i=0 ^\U\=j-i 

The improved energy estimates of Proposition 4.6 then allow us to conclude that |Ayj^| < for 
some C ~ 1. As we can easily argue in a similar way for A^ 2 ) the conclusion follows. □ 

The combination of Lemma 5.1 and Proposition 5.2 immediately yields the following result: 

Corollary 5.3. Assume that E[p(Q)^] < oo for all s > 1. Then, there exists a constant C ^ 1 
(independent ofT and of the moments of p) such that, for any k >0 and any p E [0,1], we have 

k 


A^^ — A: 
i=i 


<{Cp) 


k+l 


(5.11) 

□ 


The following lemma provides useful alternative formulas for the approximate derivatives A^’s 
(which coincide with the argument of the limit in (2.20) for po = 0), showing that they coincide with 
the arguments of the limits in (2.21) and (2.22) for po = 0. 


Lemma 5.4. Assume that E[p((5)^] < oo for all s > 1. For all T > 0, the approximate derivatives 
Alp’s, j > 1, given by (5.4), satisfy the following two equivalent formulas: 


|F|=iGcF 


|F|=iGcF 


(5.12) 

(5.13) 


where both sums X]|F|=j absolutely convergent. 


□ 


Before we turn to the proof of this lemma, let us comment on the equivalent formulas (5.4), 
(5.12) and (5.13). Formula (5.4) is the natural formula that we obtain by expanding the difference 
quotient (see proof of (5.2) and of Lemma 5.1), formula (5.12) is the easiest to use in practice (see 
e.g. Corollaries 2.4 and 2.5), while formula (5.13) is the cluster-expansion formula used by physicists. 
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Proof. We split the proof into two steps. We first prove (5.12), from which (5.13) is an easy conse¬ 
quence. 

Step 1. Proof of (5.12). 

All absolute convergence issues that we need here (for fixed T) simply follow as before from Lemma 4.3 
or similar statements (based on Lemma 4.2). For the clarity of the exposition, we discard this issue 
in the proof. Let j > 1 be fixed. Separating the cases G = 0 and G ^ 0, and noting that C'^yc’ijc 
vanishes whenever G = F, the very definition (5.4) of reads 

1^1=^' gg 

+ (-1)^'+^ ^ E [(V0T + 0 • Gf{VcP^ + 0] ■ 
l^l=i 

For any |F| = j, G 0 F, G ^ 0, by (4.2) in Lemma 4.1, satisfies 

- V • (-1)I^I+V • 

SCF\G 

Testing this equation with 69(j)T (as in the proof of (4.18) in Lemma 4.4) yields 






|F|=i Gcf 
' ' G#0 


\F\=j Gcf 
' ' G#0 


+ E E E vf 


\F\=j %% SCF\G 

+ (-1)^'+1 J] E [(V0T + 0 • Gf{VcP^ + 0] ■ 
l^l=i 

Now, for G / 0, by (4.1) in Lemma 4.1 (with H = 0), solves 

Ef 0r - V • A^VdfcPr = E • Gsycv^V^f Et- 


SCG 


Testing this equation with yields 


and therefore 


4 = - E E 




\F\=j gcf ScG 

' ' G#0 


+ E E E v5f 


\F\=j OgF scF\G 


+ (-l)^'+i J] E [{Vcfr + 0 • CF{Vcj>^ + 0] , 
l^l=i 


4 = - E E E(-i)'^'®rv4T-c'5iiG\5V5p''<A? 

|F|=jGCF5CG 


+ E E E (-1)'^'®^ V5f <))r • C^iigV^ 


.F\(GuS) ,GUS' 


\F\=j OiF sff\G 


+ E E [(V0T + 0 • GaVdp^cf^ 

\F\=jGclF 
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With the change of variables G G \ S in the first term, we observe that the first two groups of 
sums cancel, so that we are left with 

\F\=j GCF 

that is, (5.12). 

Step 2. Proof of (5.13). 

Absolute convergence issues for this part of the proof (which do not straightforwardly follow from 
Lemmas 4.3 and 4.2) will be addressed at the end of this step. Let j > 1 be fixed. Formula (5.12) 
gives 

= E E [e • GgVsP^cP^] + E E \vcPt • GgVS^^cP^] . (5.14) 


|F|=jGcF 


|F|=iGcF 


= \S^^ 




By (4.2) in Lemma 4.1 (with H = 0), 6^(pT solves 

cPt-V- AV6( (Pt=Y. • GgVsP^S, 

GcF 

whereas cpT solves 

—(pT - V • A{V(pT + 0 = 0- 

On the one hand, testing these equations with (pT and (pT respectively (as in the proof of (4.18) in 
Lemma 4.4), we obtain 


4' = - 4 E [</>T5f </>t] - E ® • AV6^ cPt] 


and therefore 


E [{V(pT + 0 • AVS^ (/>t] - E ® 

|i"l=l \F\=j 

Y^[C-AV61<Pt], 


\F\=j 


^T= E 

|F|=jGcF 

On the other hand, can be rewritten as follows: 

^T= E E (-i)''^^''^''^iE[e-c'G(v4^«+o], 

|F|=iGcF SCF\G 

which yields by the change of variables S VJ G U 


(5.15) 


(5.16) 


^T= E 

|F|=jFcF 





|F|=iFcF 


{y<pT + i) 


(5.17) 


The desired result follows from the combination of (5.14), (5.16), and (5.17). Note that the sum 
defining in (5.14) is absolutely convergent by virtue of (4.15) in Lemma 4.3, and hence the sum 
X)|F|=j (5-17) is also absolutely convergent, since its terms have just been rewritten but are still 
the same. Likewise, the sum in the right-hand side of (5.15) is absolutely convergent by Lemma 4.2 
(thus justifying the testing argument), so that the sum in (5.16) must also converge absolutely. This 
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finally proves that the sum in (5.13) is absolutely convergent too (which would not be clear 

a priori without performing this decomposition). □ 

5.2. Proof of Theorem 2.1 and Corollary 2.2. Let ^ G |,^| = 1 be fixed. It suffices to prove 
Theorem 2.1 and Corollary 2.2 for that fixed choice of What needs to be done is to pass to the limit 
T t oo in Corollary 5.3, and get rid of the additional assumption that E[/9(Q)®] < oo for all s > 1. 
For that second purpose, given a point process p, we introduce approximations for which all moments 
exist: more precisely, we shall construct hardcore approximations po of p, apply Corollary 5.3 for 
these approximations, and then pass to the limit in both the parameters T and 9. We split the proof 
into five steps. 

Step 1. Hardcore approximations of p. 

Let 0 > 0 be fixed. In this first step, we construct hardcore approximations pe of the stationary 
point process p in the following sense: for any 0 > 0, pe is an ergodic stationary point process on 
such that pg d p and pe{Q) < 9 a.s., and moreover pg p locally almost surely as 0 f oo. For 
any 0 > 0, we choose a measurable enumeration pg = (^^)n- We then define Aq as the coefficients 
obtained when replacing p by pg in . Similarly, we define (j)^g the approximate corrector and ^ 

the approximate homogenized coefficients associated with Aq ,A^q'^ instead of A^ ,A^P\ We then also 
prove the following convergence properties, which will be crucial in the next step: for fixed p G [0,1] 
and T > 0, we have 

E[|V(4^) ^0, (5.18) 

and therefore 

|4p)_4p)|^0. (5.19) 

We first give a possible construction of such an approximating sequence {pe)g. Consider the mea¬ 
surable enumeration p = {qn)n, choose independently a sequence {Un)n of i.i.d. random variables 
that are uniformly distributed on (0,1), and consider the decorated process (g„, Un)n- We then build 
an oriented graph on the points {qn, Un)n in x [0,1] as follows: we put an oriented edge from 
{q,u) to {q',u') whenever {q + ^Q) n {q' + ^Q) ^ 0 and u < u' (or u = u' and q precedes q' in the 
lexicographic order, say). We say that {q',u') is an offspring (resp. a descendant) of {q,u) if {q,u) is 
a direct ancestor (resp. an ancestor) of {q',u'), i.e. if there is an edge (resp. a directed path) from 
{q,u) to {q',u') in the oriented graph constructed above. We now construct pg as follows. Let Fi be 
the set of all roots in the oriented graph (i.e. the points of Vq without ancestor), let Gi be the set 
of points of Vo that are offsprings of points of Fi, and let Hi = Ti U Gi. Now consider the oriented 
graph induced on (q„, Un)n \ and define F 2 , G 2 , L /^2 in the same way, and so on. By construction, 
the sets F) and Gi are all disjoint and constitute a partition of the collection (q„, Un)n- Finally define 
pg := TTidJ^F)), where tti is the projection on the first factor, 7ri(q,u) = q. We easily check that pg 
defines a stationary point process on and satisfies the required properties. Ergodicity of pg easily 
follows from that of p exactly in the same way as for the random parking measure in [17, Step 4 of 
the proof of Proposition 2.1]. 

It only remains to prove the convergence property (5.18). For that purpose, we write the equation 

(p) (p) 

satisfied by the difference \ 

- cpP) - V • A<f^v{4l - = V • (4"'^ - 7l(?'))(V04 + 0- 

(p) (p) 

Testing this equation in probability with (p}fg — (pj. itself yields 

E[\Vi4l - ^T^)\^] ^ + !)]• 
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By assumption, A(0) and ^^(0) only depend on p via the restriction p\b^ for some given r > 0, so 
that the same property holds by definition for Hence, for some L > 0, 

e[|v( 4^) - (iv4^)|2 +1)]. 

Now the desired result simply follows from dominated convergence and the basic energy estimate 
E[|V(/)^^P + 1] < 1, recalling that, by definition, we have, almost surely as 0 f oo, 


Step 2. Reduction by regularization. 

In this step, we prove Theorem 2.1 and Corollary 2.2 provided we have that, for fixed T and under the 
additional assumption E[p(Q)^] < oo for all s > 1, the map p e-^ satisfies, for any po £ [0,1] 

and any A: > 1, for all —po < p < 1 — Po, |p| < l/C'po) 


1=1 


(po),j 

T 


<(pC, 


\ A;+1 


Po) 


(5.20) 


for some constant Cp^ ~pp 1, where the are equivalently given by the arguments of any of the 


(po),i| 


< for all j > I (uniformly 


limits (2.20), (2.21) and (2.22), and further satisfy the bounds |Aj, 
in r,po and the moments of p). 

Let Po G [0,1] be fixed. Consider the approximations pe introduced in Step 1, and apply (5.20) with 
p replaced by pQ (where obviously all moments of pg are finite). For any A: > 1, it follows from (5.20) 

that the map p eA ^ • A^'^g^ is smooth (on the whole of [0,1]), and a Taylor expansion of the map 
around po up to order k gives, by Lagrange’s remainder theorem, for all —po < p < 1 — po, 


k+1 


sup \Aj,g 
«e[o,i] 


ipo+up),k+l< 


<{Cp) 


k+1 


(5.21) 


« ■ Afpp'k - ( ■ AfS( - < P 

1 = 1 

From (2.7) and (5.19), we learn that 

- f ■ = i™ «■ = f ■ 474’f -«■ 474- ( 5 - 22 ) 

Hence, in order to pass to the limit T, 0 f 00 in (5.21), it is enough to prove that the limits 

^(po),l ._ (5.23) 

Ttoo 6»too 

all exist in M, for all j > 1. The combination of (5.22) and (5.23) indeed yields that for any A: > 1, 
for all —Po < p < 1 — Po, we have 


1=1 


<iCp) 


k+l 


which is equivalent to the analyticity statement of Theorem 2.1 (with convergence of the Taylor series 
at Po for all perturbations p of magnitude |p| <1/(7, —po < p < 1 —Po), and the derivatives 
are then given by the desired well-defined limits stated in Corollary 2.2. In the particular case when 
the process p has all its moments finite, the regularization in 9 can be omitted (so that only the limit 
in T remains). The proof of formula (2.23) in Corollary 2.2 is postponed to Step 5. 

We prove (5.23) by induction. The proof of the statement for j = 1 is similar to the proof of the 
induction step, and we only display the latter. Assume that the limits A^^AJ = lim^ lime 
exist in M for all 1 < j < fc, for some A; > 1. We shall then prove that the limit = 

lim(p limp also exists in M. As ^ is bounded uniformly in T,9, it converges to some 

limit 4^“^ G M as 0 t 00 up to extraction. Passing to the limit 0 f 00 along a subsequence in 
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inequality (5.21) with k replaced by /c + 1, and using the induction assumptions and (5.22), we obtain 
for any -po < p < 1 - po, 


e • - C • A 


(po) 


^ p^ lim A, 
i=i 


(po),j 

T,6» 




<{Cp) 


k+2 


This proves that L^p°') satishes 


k 

4''"'= lim (p-*-'(£.4«+'’>{-O/l^>0-Ep-'-'limA«-4 

-Po<p<i-po ^ i=i ^ 

where in particular the limit must exist. Since the right-hand side does not depend on the extraction, 
is uniquely dehned, and = lim^ does exists in R. A similar argument for the limit 

in T shows that limy lim^ exists in R, so that (5.23) is proved. 

Step 3. Reduction by restriction to po = 0- 

Let T > 0 be hxed and assume that E[p(Q)^] < oo for all s > 1. In the present step, we prove 
that it suffices to check the result (5.20) at po = 0: more precisely, it suffices to show that the map 
p I—^ satishes, for all /c > 1 and p G [0,1], 


j=i 


<{Cp) 


k+l 


(5.24) 


for some constant C ~ 1, where, for any j > 1, A^ is equivalently given by formulas (5.4), (5.12) 
and (5.13), and satishes the bound |Ay| < for all j > 1, uniformly in T and the moments of p. 

First consider po G [0,1) and positive perturbations po +p with p > 0. For 0 < p < 1 —po, choose a 
sequence {dn°'^^)n of i.i.d. Bernoulli random variables (independent of all the others) with parameter 
= 1] = p/(l — Pq)i and consider the twice perturbed coefficients 

^(Po,p) = ^ + (1 - 

n 

= + E (<1 - ‘'(J*'”’)*' - Ij,. 

n 

The held A^p°’P'^ has by dehnition the same distribution as A^p°~^p\ and it is a perturbation of A^P°^ 
with perturbation parameter p/(l — po) (and with perturbed medium A'). Applying to A^p°’P^ the 
result (5.24) around 0 (which is assumed to hold), we deduce that the map p i—?> 
satishes, for any A; > 1, for all 0 < p < 1 — po. 


^- ? • a: 


(po) 


«-E 


i=i 


_A 

{1-Poy ^ 


(po)j 


< 


Cp \ 

1 -PO/ 


fc+1 


(5.25) 


where, for any j > 1, is the j-th right-derivative at 0 of the map p i—)■ corresponding to 

the “reference” coefficients A := A^P°'^ and the “perturbed” coefficients A'. The cluster formula (5.13) 
reads in that case 


^4o)j 


\F\=j Lgcf -I 


(5.26) 


where the sum is absolutely convergent. Now note that the argument of the expectation vanishes 
whenever F n ^ 0, while, otherwise, if T n = 0, the argument of the expectation equals 

^ (-l)inGI^ . ^GU(E(ro)\F)(y^GU(E(Po)\F) ^ 

G(ZF 
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As this expression is obviously independent of the event [F n = 0], we can then rewrite 


\F\=j 


^FnE(po)=0 ^ • A 

GCF 


E(^o)uG(v<^E(po)uG^^) 


P[F n = 0 ]E Y (-1)'^'^^'? • + 0 


\F\=j 


'-GcF 


= ii-PoyAi^ 


{Po),k 


(5.27) 


where is defined as the argument of the limit (2.22). The expansion (5.25) then becomes, for 

any A: > 1, for all 0 < p < 1 — pq, 


i=i 


< 


Cp \ 
1 -Po/ 


fc+i 


(5.28) 


Moreover, recalling the bound |A^| < for the right-derivatives at 0, which is assumed to hold 
for any choice of the coefficient (the constant C only depends on R,V,d,\), we conclude, for all 
Po e [0,1), 

< C^(l -po)"^- (5.29) 

Note that this estimate for the derivatives deteriorates when po gets closer to 1. This difficulty is 
overcome by considering negative perturbations, that is, looking at left-derivatives, as we do now. 

Let us now consider po G (0,1] and negative perturbations at that point. For 0 < p < po, 
choose a sequence ((i^°’ ^'^)n of i.i.d. Bernoulli random variables (independent of all the others) with 
= 1] = p/poj and consider the twice perturbed coefficients 

^(P0,-P) = + Y {d^t'-^^A + (1 - 1 

n 

= j. + E + (1 - - hFyF +- diy’-FA') 1 j.. 

n 

The field has by definition the same distribution as and it is a perturbation of 

with perturbation parameter p/po (and with “perturbed” medium A, instead of A'). Applying to 
^(po,p) result (5.24) around 0 (which is assumed to hold), we deduce that the map p i—)• = 

^ ■ A^° satisfies, for any A: > 1, for all 0 < p < po. 


« . Ap-'k - ( 


Apk-T. 

J = 1 


Po 


< 


Po J 


k+l 


(5.30) 


where, for any j > 1, A^°'^d jg j-th right-derivative of the map p i— A^\ corresponding to the 
“reference” coefficients A := A^'pA and the “perturbed” coefficients A. The cluster formula (5.13) gives 
in this case 




\F\=j ^GCF 


where the sum is absolutely convergent. Arguing as above, the argument of the expectation vanishes 
unless F C hence, by the independence assumption, we obtain 


^{P0),j 


\F\=j ^GCF 
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or equivalently, by the change of variables F \ G H, 



(-i)Vo E ® 

l^l=i 


- HCF 




where, as before, is defined as the argument of the same limit (2.22). The expansion (5.25) 

then becomes, for any A: > 1, for all 0 < p < po; 


^^. 4°)^ _ ^ 

f=i 


< 



(5.31) 


By the bounds \AY\ < for the right-derivatives at 0, we conclude, for all pq G [0,1), 

< (74^'. (5.32) 

Combining (5.28) and (5.31) then directly yields the desired result (5.20). Moreover, combin¬ 
ing (5.29) and (5.32) gives, for any J > 1, the uniform bound 

|^^o)J| ^ mm{C^P qY — Po)~^} < (2(7)-^. 

Finally, arguing as in Lemma 5.4 (where the argument is performed at po = 0 and proves the 
equivalence between formulas (5.4), (5.12) and (5.13)), we see that, for fixed T,6, the cluster formula 
for A'Y°'^’^ , that is the argument of the limit (2.22), is equivalent to the formulas given by the argument 
of the limits (2.20) and (2.21). 


Step 4. Conclusion. 

Let T > 0 be fixed, and assume that ]E[/9(Q)^] < oo for all s > 1. For any A: > 1, Corollary 5.3 
exactly asserts (5.24), Lemma 5.4 ensures that the A^’s are equivalently given by formulas (5.4), 
(5.12) and (5.13), and Proposition 5.2 gives the uniform bounds |A^| < , for all j > 1. By the 

previous steps, this proves Theorem 2.1 and Corollary 2.2. 


Step 5. Exact formula for the first derivative. 

In this last step, we further assume that p{Q) < Bq a.s. (so that in particular all the moments 
are bounded, and we can thus everywhere omit the regularization in 9), and we prove under that 
assumption the validity of formula (2.23) in Corollary 2.2. More precisely, we need to prove that we 
can pass to the limit in T inside the formula for the first approximate derivative 


Ab>o),i 






i.e. we prove that the well-defined limit A^^°)7 


lim^ Ay 


is given by the following formula: 




(5.33) 


with an abolutely converging sum. As before, we can restrict to pQ = 0, and shall prove that the limit 
A^ := limy Ay exists and is given by 


A^ = ^E 


(V((. + o-c^''^(v4^ + 0 


We start by showing that the sum is absolutely convergent. Decomposing -|- V(j), 

using assumption (2.1) in the form < 1 (see also (4.20)), and recalling the elementary energy 

estimate E[1 -|- |V(/)p] < 1 (see (2.5)), we have 

|Ai| < ^E[lj„(l + \VcPY + |V4^12)] < 1 + ^E[|V,5^"4|2], 
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where the last sum is hnite by Lemma 4.4 (for fc = 1) and the fact that weakly 

in LL(M'';L 2 (L!)). 

We now prove that limj’ = A^. Given L ~ 1, the additional assumption p{Q) < 6o a.s. implies 
by stationarity p{Bji^L) < COq =: Z, with C ~ 1. Hence, we can choose the measurable enumeration 
{Qn)n of the point process p in such a way that n {qn)n C {qn)n=i- Defining 


al :=E 


-JBl 


and 0 ^ 2 . ■= ® 


JBl 


we observe A^ = aj and A^ = 'Yln=i^TL- Indeed, by stationarity (together with absolute 

convergence), A^ = aj, so that A^ = 'Ylin=i by the choice of the measurable enumeration, 

and likewise for A^p. Therefore, it is enough to prove limj^a^^ = for any 1 < n < T. Since 

V(j)ip^ —^ V())I”I weakly and V(j)T —t Vi/) strongly in Lj^q^(]R‘^,L^( n)) (see [11, Theorem 1]), we 
directly get ^ aj as T f oo, for any n, as desired. 


5.3. Proof of Corollaries 2.3, 2.4 and 2.5: Clausius-Mossotti formulas. In this section we 
further assume that K[p{Q)‘^\ < oo. (Note that, in the case of Corollaries 2.4 and 2.5, this directly 
follows from assumption ( 2 . 1 ) together with the fact that we are then dealing with ball inclusions of 
fixed radius.) 


5.3.1. First-order universality principle. Set — |J^gg{p) Jn- The volume fraction Vp of the per¬ 
turbation is defined as follows: 

E[|LQn J(D|] 


v„ := lim 
^ Ltoo 




or equivalently, by stationarity of the inclusion process, 

Up=limL-'^ ^ E[|(z + Q)n =E[|Qn 

L'\oo 

zGLQDZ^ 

An inclusion-exclusion argument gives 

^[^nGE(p'>\Q D Jn\] ~ ®[I^n,meE(p) IQ G H Jm\] ^ Vp A ®^[I^neE(p) IQ ^ Jn\]- 

n n^m n 

By the independence between the Bernoulli process and all the other random variables, this 
turns into 

P^E[|Qn Jn|] ^ E[|(5 n Jn n Jm|] < Vp < P^E[|(5 n Jn|]. 

n n^m n 

As Jn C Bji{qn) for all n, we note that, by the assumption K[p{Q)'^] < oo, 

n Jn n Jm|] < E ^qn,qm&Q+Bji = '^[p{Q + Bpi) ] < OO, 

n^m ^ n^m 

so that we have indeed proven 

Vp = P^mQ n Jn\] + 0{p^) =:p-f + 0{p‘^). 

n 

If 7 = 0, then UnJn = 0 so that = Ahom and Vp = 0, and there is nothing left to prove. If 

7 / 0 , since 


7 < E 






= E[p{Q -\- Br)] < 00 , 


we have Vp ~.y p -|- O(p^). In particular, the expansion in p in Theorem 2.1 at first order can as well 
be rewritten as an expansion in Vp'. at hrst order at pQ = 0 , we have, for any p > 0 , 

— A. _L JO'jI _L O) f'I^ 


,2 
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where K is given by 


e-i^e = -e-4ome = -]E 

^ llCUll ^ 




If in addition the random volumes | J^|’s are i.i.d. and independent of the point process p (and of 
its enumeration), then 7 can be computed more explicitly. By stationarity of the inclusion process, 
for any L > 0 , 


7 = lim L 
Ltoo 


^\LQr^Jn 


where we can estimate 


E 


E \LQ n Jn\ < Y,n\..^LQ\Jn\] = E[| Jo° |]E[p(LQ)], 


and also 


E 


E \LQ n Jn\ > Y.^[t,„^iL-R)Q\Jn\] = E[|Jo°|]E[p((L - R)Q)]. 


(5.34) 


(5.35) 


Now, for all continuous and integrable functions / : —)■ R, we have E[^^/(g„)] = E[J’ fdp] = 

f fdK[p]. Since p is stationary, the Borel measure E[p] is translation-invariant, and hence, since it is 
locally finite by dehnition, it is a multiple of the Lebesgue measure: E[p] = adx for some constant 
a G R"*", which is characterized e.g. by cr = E[/9((5)]. In these terms, (5.34) and (5.35) give 


aE[| Jo°|] = E[| Jo°|] limL-'^E[p((L - R)Q)] < 7 < E[| Jo°|] lim L-'^E[p(LQ)] = aE[\JS\], 

1j Jj 


which means 7 = (tE[| Jq |], and thus 

The matrix K then takes the form 

i-Ki = 


Vp =pcjE[|Jo|]. 


(5.36) 


1 


aE[|J°|] 


E 




Further assuming that p is independent of A, of {A'^)n (as well as of the random volumes | J°|’s), we 
note that the random variable E[(V(/>(0) -|-.^) • C'’f”^(0)(V(/)’f"’^(0) +0llp] o^ly depends on the point 
process p through the point so that it can be written as f{qn) for some measurable function /. In 
these terms, we get 


i-Ki = 


aE[\J°\] 


E 




E[J°|] 


[ f{x)dx, 


which does clearly no longer depend on the choice of the point process p. This proves Corollary 2.3. □ 


5.3.2. Electric Clausius-Mossotti formula. We consider the case when the inclusions are spherical 
Jn = Bn{qn), and the unperturbed and perturbed coefficients have the form ^ = aid and A' = (did 
respectively. We shall compute explicitly the hrst derivative f of th® perturbed homogenized 

coefficient at 0, as given by formula (2.23). As inclusions are balls of fixed radius R, assumption (2.1) 
implies p{Q) < 9q a.s. for some constant 9q > 0, so that we can indeed apply formula (2.23). 

Since A is constant, the unique gradient solution of —V • A{S/(j)^ -|- ^) = 0 is clearly = 0. Let 

now n be fixed. The solution G Ff^(R'^) of —V • + ^) = 0 is easily checked to be 

unique if it exists, and its existence follows from a direct computation. Since is constant both 
inside and outside Bji{qn), the solution is radial and of the form 


='4)^{x-qn) 


C{x-qn)-f, ioi\x-qn\<R] 

Qt lx-qnH 1^ _ I 
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SO that its gradient satisfies 


-qn) = i g/ / 


for lx — qn\ < R', 

_ x-g \ I _ I 

\x-qn\ \x-qn\J ^ I 


Since is radial and in it is continuous, which implies that C' = The normal 

component of must also be continuous through the sphere, so that we conclude 

a — d 


C = 


(5.37) 


j3 + a((i — 1) 

This allows us to turn (2.23) into an explicit formula for the first derivative : 

? • = E + 0] 

n 

= ]E[^ (/3 - a)lB«(,.) (Ce + 0] = (1 + C)(/3 - a) E [ ^ 


From Paragraph 5.3.1 above, we learn that E[^^ /(^n)] = ^ f f{x)dx for any continous and integrable 
function /. Hence, 


^4ome = (l + Cl)(/3-a)E 

so that expression (5.37) for C yields 




= a\BR\{l + C){/3 - a), 




ad{(3 — a) 
(3 + a{d — 1) 


In the present case, formula (5.36) holds true and gives Vp = pa\BR\. The conclusion of Corollary 2.4 
now follows from Theorem 2.1. □ 


5.3.3. Elastic Clausius-Mossotti formula. We consider the case of spherical inclusions J„ = Bii[qn) 
and assume that both the unperturbed stiffness tensor A and the perturbed stiffness tensor A' are 
constant and isotropic — we denote by K, G and K', G' their respective bulk and shear moduli. We 
shall compute explicitly in that case the first derivative ^ of the perturbed homogenized 

stiffness tensor 4E at 0, as given by formula (2.23). Indeed, as inclusions are balls of fixed radius R, 
assumption (2.1) implies p{Q) < Bq a.s. for some constant 6q > 0, so that we can apply formula (2.23). 

Let G be symmetric. Since A is constant, the unique gradient solution of —V-H : = 

0 is clearly = 0. Let now n be fixed. As shown e.g. in Section 17.2.1 of [32], equation 

— V • + f,) admits a (necessarily unique) solution in iL^(M'^). Inside the inclusion 

RniQu), that is, for all |x — qn\ < R (see equation (17.84) of [32]), 


V0f^(x)+C = ld 


TifK + d 

d K' + /3 


+ 


e-Id 



G + a 
G' + a 


(5.38) 


where a, d are defined by (2.25). Recalling that : A : x = : y + ATr^Try for any symmetric 

y G we can now explicitly compute formula (2.23) for the first derivative 


It . 4(0)4 . e — 1 

2^ * "^hom * ^ 2 


wl”'+o 


E ^ "B-h.) ((O' - G)« : (vd”' + 0 + 1(A' - A) Tr4 Tr(V,^f > + {) 
















ANALYTICITY OF HOMOGENIZED COEFFICIENTS UNDER BERNOULLI PERTURBATIONS 


45 


and hence, using (5.38), and recalling from Paragraph 5.3.1 that 
1 


2cr\BR\ 


!.1‘ : { = ((G' - G) + - >~)) 

+ (|{|"-l(rV{)=') (G'-G) 


G (X 
G' + a' 

In terms of bulk moduli K = X + 2Gjd and K' = X' + 2G' jd, this takes the form 


_ C ■ /li^)d . c — tl2/' 


2a\Bii\ 




K + ^ 
K' + P 


1 


+ (lCr-^(TrO^)(G'-G) 


G ot 
'G' + a 


(5.39) 


As formula (5.36) again holds true in the present case and gives Vp = pa\Bii\, Corollary 2.5 now 
follows by Theorem 2.1. □ 

5.4. Proof of Corollary 2.6: convergence rates. Let po £ [0,1] be fixed, and assume that 
E[p(Q)®] < oo for all s > 1. Estimate (5.21) in the proof of Theorem 2.1 (see Section 5.2) then yields, 
for all A; > 1, 


(po),i _ aIpoI’T 
2T ^ 


i=i 


< (Cp)^+1 + - A^^^l 


< (Cp)"+1 +E[|V(4^°+^) - +E[|V(4^°) - 4??))|], 

and hence, combining this with assumption (2.26), 


1=1 


ipo),j _ A (po)jl 
T ^2T V 


< (Cp)^+i + C7(r). 


By induction, we easily see that this implies, for all J > 1, 

< {2Gy+^-f{Tf~\ 


(5.40) 


(5.41) 


Estimate (5.40) with A: = 1 gives < G'^p + G'y{T)/p, which turns into (5.41) for 

j = 1 with the choice p = 7 (T) 2 . Assume now that the result (5.41) is proven for all 0 < j < J. 
Then, equation (5.40) for A: = J + 1 gives 


l^bo),J+i _ < (jJ+2p ^ ^ ^(2C) V^"S(^) 

1=1 


2-J 


With the choice p = 7 (T)^ , and noting that (I + 1)2 < 1 for any I £ N, this turns into 


- ^ +G-i{T) 

J 

1=1 


l-(J+l)2--^-l 
2-l(l-(J+l-j)2-(-^+l-l)) 


J 


< 


GJ+^^{Tf + G-f{Tf + ^(2C)^'+S(r) 

1=1 

^2 + ^2^+^j < (2C)^+%(r) 

which proves (5.41) by induction, and concludes the proof of Corollary 2.6. 


2-1-1 


□ 
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